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Introduction. 

Lame's functions or ellipsoidal harmonics have been successfully used in many 
investigations, but the form m which they have been presented has always been 
such as to render numerical calculation so difficult as to be practically Impossible. 
The object of the present Investigation Is to remove this imperfection in the method. 
I believe that I have now reduced these functions to such a form that numerical 
results will be accessible, although bv the nature of the case the arithmetic will 
necessarily remain tedious. 

Throughout my work on ellipsoidal harmonics I have enjoyed the immense 
advantage of frequent discussions with Mr. E. W. Hobson, He has helped me 
freely from his great store of knowledge, and beginning, as I. did, in almost complete 
Ignorance of the subject, I could hardly have brought my attempt to a successful 
issue without his advice. In many cases the. help derived from him has been of 
Immense value, even where it is not possible to indicate a specific point as due to him. 
In other cases he has put me In the way of giving succinct proofs of propositions 
which I had only proved by clumsy and tedious methods, or where I merely felt sure 
of the truth of a result without rigorous proof. In particular, I should have been 
quite unable to carry out the Investigation of § ! 9, unless he had shown me how the 
needed series were to be determined. 

My original object in attacking this problem was the hope of being thereby enabled 
to obtain exact numerical results with respect to M. Poincaiie's pear-shaped figure 
of equilibrium of a mass of liquid In rotation. # But I soon found that a partial Inves- 
tigation with one particular point in view was Impracticable, and I was thus led on 
little by little to cover the whole field, in as far as it was necessary to" do so for the 
purpose of practical application. This paper has then grown to such considerable 
dimensions that it seemed best that it should stand by itself, and that the discussion 
of the specific problem should be deferred. 

A paper of this kind is hardly read even by the mathematician, unless he happens 
to be working at a cognate subject. It appears therefore to be useful to present a 
summary, which shall render it possible for the mathematical reader to understand the 
nature of the method and results, without having to pick it out from a long and 
complex train of analysis. Such a summary is given in Part 111. 

PART I. 

Formation o:ir the Functions. 

§ 1. The Principles of Ellipsoidal Harmonic Analysis. 

The basis of this method of analysis is expounded In various works on the subject, 
begin with a statement of results in my own notation. 

* A paper giving the required result will be presented to the Society in the autumn, —\July 2, 1901.] 
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If ^A u£, Uo 2 denote the three roots of the cubic 



i ? ov z -> 



o 



?/' 



S 3 



0,0 7 o . n | o , o 



W -\ 



•f U : 



it may be proved that 



ar 



(a 3 4- ?0 {a? + V) (a 2 4- n*) 

(& 3 — a 2 ) (c 3 — a 3 ) 



and ?/ 2 , a; 2 may be written down by cyclical changes. 
If for brevity we write 

A,, 3 = n* + a\ B/ = u,? + b\ C,~ - »/ + c 3 , (ft - 1, 2, 3), 



Laplace's equation becomes 






rf \« 



d 



\9, 



<) ( A ^ B A ;,k ) v < + ("» 9 - «i 3 ) ( ^ c » ^ ] v < 



(V f 5 .' 



+ «-<-) A 3 B 3 C 



r7 '^ 



The solution is 



V,- = TJ.XLTL, 



where U l5 IJ 25 U 3 are functions of u ly t%, % respectively, and satisfy 



rl \s 



~i{i + iW + « 8 ]U X , 



and two other equations with suffixes 2 and 3, involving the same /c, a constant, and 
the same 2, a positive integer. 

If a, 6, c are in ascending order of magnitude we may suppose u L 2 to lie between 
— a 3 and oo , u£ between — c? and. ~ (V 2 , and % 2 between — ?/ and —a 2 . 

If s l9 5 2 ., 6 f 3 denote the three orthogonal arcs formed by the intersections of the 
three orthogonal quadrics, 



( ds { 
u \du\ t 



2 



( u £ - y) «- - U*) 



and two other equations found by cyclical changes of suffixes. 



§ 2. Notation; limits of /3 so as to represent all Ellipsoids. 
I now change the notation, and let the three roots be defined thus : — 



o 
US 



7 9 9 



o 7 9 

Ct'O — IV 



7 o o 

1 — j3 cos 2$ 



o 



1 - /3 



where ^ ranges from oo to 0, jx between f 1, ^ between and 27r. 
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Let the axes of the fundamental ellipsoid of reference be 



1/ , «- / t , ? 

c 2 = 0, 

The ellipsoid defined by v has its three axes a, b 9 c given by 

/ 1 + /3\ 
a 2 = & 2 f p 2 — t- r) , 6 2 = F(i/ 3 — 1), e 2 = &V, (a < & < c). 

This mode of defining the axis is such as to indicate the relationship to the 
prolate ellipsoid a = 6 < c. But another hypothesis may be made which will brine 
the axes into relationship with those of the oblate ellipsoid e& = & > c ; for if we 
take a new h, numerically equal to the old one but imaginary, and replace v 2 by — £ 2 , 
we hare 

a 2 = ¥ (? + y~;) , A 2 = P (£ 2 + 1), c 2 = F£*, (a > b > c). 

If /3 be made to range from to oo , all possible ellipsoids are comprised in either 
of these types. It will, however, now be shown that, by a proper choice of type, all 
ellipsoids may be included with the range of /3 from to |. 

Let us suppose the axes to be expressed in three forms, as follows : — - 

(1.) ■ (2.) (3.) 



W = **(„* - 1) = V(£i s + 1) = * 2 *L* 






e 2 



Then we have 



Therefore 



hV = b»(c* + L±-|) = vw - i). 



6 * — 1 — /3 — ^ "~^ 3 1~ A 

O 7 7 /-,7ti AJi i-//U Dg 

q* —«. (^ ^^ ^* ^^ - i -.:.-. x - .-^ 

£ 3 - a 2 2/3 1. -ft 1 -f- /3 3 



o 



- & 2 1 - /3 2/^ 2/5 3 J 

c 3 — a 3 1 4- j3 1 + /^ 1 — ft? 



O 7 O 



£/ X ~~" /.!> Zifo-i *jkjq 



2c» - a 3 - 6 3 P ~~ 1 + ; ;, »A ~ 1 '- 3/3 



PROFESSOR G. H. DARWIN OX ELLIPSOIDAL HARMONIC ANALYSIS. 465 

Now let /3 increase from to c© . 

As /3 passes from to |- ? form (1) is appropriate. 

As yS passes from ^ to 1, /3 t decreases from |- to 0, so that form (2) is appropriate. 

Lastly, as /3 passes from I to co , /3 2 increases from to ■§-, so that form (3) is 
appropriate. 

But we might equally well have written forms (l) and (3) so as to involve £, and 
form (2) so as to involve v, and it follows that all possible ellipsoids are comprised in 
the range of /3 from to \, provided that the type be appropriately chosen. 

The developments in this paper are made in powers of /3. It will, therefore, be 
well to show that there is a class of ellipsoids, analogous to ellipsoids of revolution, 
which might form the basis of developments similar to those carried out below. 

Ellipsoids of revolution are defined by the condition 

a 2 — c 2 = b 2 — c 2 , or a 2 = & 2 . 
In the class to which I refer 

a~ — <r ■ = er — 6% or <r = -&■ (a- + '> )• 

Ellipsoids of this kind are given by ft = /3, = — /?., = -J- ; for in tin's case 
lr = Ucr + ^ 2 ). Thev are also given bv 

ft = ~o 9 - p } = fr 2 = i. f or then c 2 = |- (a 2 + 6 2 ). 

Hence if we only allow /3 to range from to £, /? = corresponds with ellipsoids 
of revolution, to which spheroidal harmonic analysis is applicable ; and /3 = £ 
corresponds with this new class for which the corresponding analysis has not yet 
been worked out. 

"We shall see below that the solid harmonic for this case where ft = % will be of 
the form B (V) B (/x) E (<£), where B and E satisfy the equations 

( r 3 + j) ^ _ -, } U + 2j , 3 | ] _ -^ + ^ ^ B + ^ B = (h 
cos 2<f> -~ - sin 26 ~ + { ({ -f 1) E cos 2<f> - srE = 0. 

I am not clear whether or not it would be advisable to proceed ah initio from these 
equations, but at any rate I shall show hereafter how the B- and E-functions may be 
determined from the analysis of the present paper with any degree of accuracy 
desirable. 

If it were proposed to use the functions corresponding to /3 = -J as a basis for the 
development of general ellipsoidal harmonics, we should have to assume 
yol. cxcvu. — a. 3 o 
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a 2 = Jcf*i/' 2 9 W = k hl (v f% — 1), c 3 = k' % [v'*— - / , 



2_ 



or else a 2 = & /2 £ /2 ? W = V* (£'* + 1), c 2 = H^ + ;"- - 

\ 1 — 97/ 

The developments would then proceed by powers of 77. 

In order to discover what is the greatest value of tj which must be used so as to 
comprise all ellipsoids, when we proceed from both bases of development, a com- 
parison must be made between this assumption and the previous one. Suppose in 
fact that 

a 2 = l$h> - ™j|) = *'*£'*; h* = P(i^ - 1) = &'*(£'* + 1) ; 
Then W - a 2 = ~ -- = jfc* ; c 8 - fe 2 = fc* = F 3 i-^ 

1 — - P 1 — 1/ 

1 , , p 2/3 1 — 97 1 ~~ 3/3 

ana therefore : — -— = , or 77 ~ - — -- , 

1 — yS 1 + 97 ' 1 + /3 



When 77 and /3 are both equally great, they must each equal the positive root of 
1 - 3/3 

1 + ft 

1 1 



/3 = - — ---,-. This root is A /5 — 2 or *236, Thus the greatest values will be 



J3= v 



5 + 2 ' 4-236 



In this case rf = /3 3 — -^ very nearly, whereas when ft = |-, /3 2 = -3-. Thus if the 
developments were to stop with ft % we should double the accuracy of the result. 
However, I do not at present propose to carry out the process suggested. 

§ 3. The Differential -Equations, 

X -~ /3 cos 2(f) 
We now put u^ = &V, u£ = P/x, 3 , w 3 3 = F . -- — - ;~-~ ; 

.1. "™~~~ p 

a a _ _ 1% L±l ft s _ _ w c a _ . 



1 -£ 
and find from the formulae of § 1 , 

7« 3 i + p v 

?/ 3 



J. "T* 


-w- 


1 


4- 


e\ 


cos 


2 <£> ! 


1 - 


I 


__ 


p) 


- 1 


) sin 2 (j) , 













(v % — 1) (/x 2 — l)sin 2 <£, )- . . . (1). 



£ 2 „ ; , 1 — /3 cos 2(f) 



* 3 r 1 + /3 j 
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It will be observed that y is independent of /3, and that it has the same form as in. 
spheroidal harmonic analysis when /3 vanishes. Since /x s is less than 1 and v 1 greater 

than r-—-z , x and y are real. 
1 - ft J 

In all the earlier portion of this paper I always write jut — 1 and not 1 — /x 2 , so 
as to maintain perfect symmetry with respect to v and \l. 
We now have 



A 3 _ 
•^1 — 


k' z 


('" 


■'■ + /S'\ 

T-p) ' 


B x * 


= tf(f 


.— 


1). 


(V- 


7 o 


A3 — 

xjlo — — 

(4/ 


& 


(m 3 - 


1 + /3\ 
~ 1 - p) ' 




= F'(/A 3 


— 


!)» 


c/ 


= P/x, 3 ; 


A 3 _ 
ii. 3 — 


— - 


2F/3 
1 - 


cos 3 $> 
' 




2W/3 sin 3 
~ 1 - P 




c 3 3 


, o 1 — /3 cos 2<£ 


nst dfvr 


int 


A ttlP 


' rliflfiftrftntia.l < 


"vnfimtr 


■*rfi invnli 


7W 


1 in nn 


r fimifl. 


f.inns f,Ln« • — - 



D- 



(1 -fi) 



x AjBjGi a 



Then 



Al^i^i 



<Z 



/^ x c^ * 



D 



(1 



/J) x 



t^^j (1 — /3) 



— rttf ? 



D, = (1 - m^^A 



2 



Do 



y-l.(l - pf 



x iio-DqOq Lv 



jS)*(^~ j-z-|)V 



D 2 = (1 - £)*(/ti 2 

\ 

\ 

D 3 = (1 — /3cos2<£) 



i 



xjLoXJo v-/o 



i + fl y, „ 

d 



1) 

1)' 



dfi* 



•\ 



y 



/i/JL/c) 



^%f/w 3 (1 — /3) s 



d 6 6 u z du z 



\ 



Hence our differential equations are 



- IX 3 IL 

o_ 6 

1 -£~ 



»(»+ J> 2 +'^ 



# 



Ji^fS 

/vivo Lotto 



# * • e « 






H)< 



(1 - 0? 



U 1? a similar equation with suffix 2, and 



i (i + 1) 



1 



fiG0S2(j) K" 



P 



k 






u 3 . 



(2). 



Let us replace «r 3 bv another constant such that 

(i(i + l)i/ 3 + -~)(1 - £) = *>* + 1)[> 5 (1 - J3) ~~ 1] + .v 2 - /3a- , 



so 



AC" 

F 3 






3 O 
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In this formula s is a constant integer and or a constant to be deterniined. 



Our equations are now 



) 



(i + 1) |^(i — /8) — I ] — .r -f- /3crJ L 



J \ 



o, 



-A 



a, similar equation for /x 



r 






and 



I)/ — i(i + 0/3cos2</> + s* — /3o-'\XJ, 







And Laplace's equation is 



fxr 



-£f) w + ( : 



1 



/3coa 2^ 



^ 



I)^ _ („2 _ ^a) D 3 a 



ILILIL = o. 

O. .V O 

(4). 



Laplace's operator V 2 is equal to the differential operator in (4), divided by 



¥{/ 



ox , 1 - /3cos2c/)\ /l - /3cos 2cf> 
/x~) ( i/- 



F 



i - # / v 1-/3 

It is well known that in spheroidal harmonic analysis there are two kinds of 
functions of v and /x which satisfy the differential equation, and they are usually 
denoted P/, Q/'. The Q-functions of the variable /x have no significance, so that 
virtually there are P- and Q-functions of v, but only P- functions of /x. The like is 
true in the present case, however, with the additional complication that each of the 
functions may assume one of two alternative forms. I adopt a parallel notation and 
write for V\ and XL either |3/, ®/, or P/, Q/, as the case may be. Since v and /x 
enter in the first two equations in exactly the same way, we need only consider one 
of them, and we may usually write simply (for example) *|jj/ where the full notation 
would be :P/(V or /x). In the early part of the investigation I shall only refer to the 
P-functions, and the Q-functions will be considered later. 

In spheroidal harmonic analysis the third function is a cosine or sine of s<j>. So 
here also we find functions of two kinds associated with cosines and sines, which I 
shall denote (£/, <£/, Oi\ S/ ? the variable </> being understood. 

Throughout the greater part of this paper the functions will be of degree denoted 
by i 9 and it seems useless to print the subscript i hundreds of times. I shall accord- 
ingly drop the subscript i except where it shall he necessary or advisable to retain it ; 
for example, |j3 s> will be the abridged notation for !$i{y). 

The operators involved in the differential equations (8) will occur so frequently 
that an abridged notation seems justifiable. I therefore write 









D 



2 



i(i+ i)|y(i -fi)-i 

i (i + 1) /3 cos 2<j> + s 



S 2 + /3(T 






/8cr 



where 



D^l-^W** 



I) 3 = (1 — yQcos 2<£)* 



I + yS \* , o ., VL d 



i 



. (5 



a 



'4 



PROFESSOR G. H. DAK WIN ON ELLIPSOIDAL HARMONIC ANALYSIS. 469 

The equations are then 

Xs(€ s or $>< or G ? or S*) = J ^' 

§ 4. 2%e Forms of the Functions* 

It is well known that the function U; is a linear function of k l of degree i made up 
in one of the eight following ways : — • 

1. When i is even, a linear function of Uf of degree \i. 

2, 3, 4. When i is odd, a linear function of u^ of degree %(i -— 1), multiplied by 

xjLi j Ol J3i , OjL Vy j • 

5, 6, 7. When i is even, a linear function of n L s of degree |(^ — 2), multiplied by 

B 1 C 1> or C 1 A 1 , or AjB^ 
8. When i is odd, a linear f auction of u^ of degree h(i — 3), multiplied by 

These eight classes may be conveniently specified by the initials O, A, B, C, BC, 
CA, AB, ABC, but it is better to rearrange them according as they are associated 
with the evenness or oddness of i and s 9 and with the cosine or sine functions. This 
new grouping may be defined by a shorthand notation involving the initials E, 
and or S, which shall denote successively the evenness or oddness of i and s, and 
cosine or sine. 

We shall see below that this arrangement is as follows : — • 

O or EEC ; i even, s even, cosine. 

AB or EES ; i even, s even, sine. 

A or OOC ; i odd, s odd, cosine. 

B or OOS ; i odd, s odd, sine. 

or OEC; i odd, s even, cosine. 

ABC or OES ; i odd, s even, sine. 

CA or EOC; i even, 6* odd, cosine. 

(IB or EOS ; i even, s odd, sine. 

Since the several functions are linear in uf, they are in the new notation functions 
of v 9j or /x 3 , or of v 9j — 1 and jm z — 1. 

Hence ffi(v) and P s (v) involve linear functions of v % — 1 of various degrees 
multiplied by various factors ; and the same is true of the functions of p. 

In the case of the third root the linear function of powers of cos 2<f> may be 
replaced by a series of cosines of even multiples of <f>. Further, in forming the C, <S, 
C 5 S functions we may regard A 3 as being cos (/>, B a as sin <p, and C ; > as (1 — /3 cos 2<£)% 
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since this only amounts to dropping constant factors which may be deemed to he 
included in the, as yet, undetermined coefficients of the several series. 

I will now consider in detail the forms of the several P-f unctions of v (those for \l 
following by symmetry), and at the same time indicate more precisely the nature of 
the notation adopted. 

In the following series, indicated by 2, the variable t is supposed to proceed from the 
lower to the upper limit by 2 at a time. The reader will be able to perceive the 
manner of the formation of the functions when he bears in mind that 



. xjL I — — /I 



'£ V 



1 - & 



k(v~ — t)% Cj — lev. 



Type O or EEC ; $' = %a t (v 2 ~ if 



Type AB or EES ; P' = 2 a* (^ - 1 )* <* - J > v % 



Type A. or OOC ; P' = 2 a t (^ - 1 ) * <' - r > ( 



i (l - 1) / j,3 



1 + j8 \* 

1 - B, 

1 +jS' 
1-8 



Type B or OOS ; $ s = 2 a, (y 3 - 1)* . 



Type or OEC 



] 



I) 



(*~ I) 



rn 



Type ABC or OES ; 



Type CA or EOO ; 



ASH^ vA* 



^ (V 2 



■|U(*-2)/„2 



^ 






I + yS 

1 +J3\* 
1 - A 



Type OB or EOS ; f s = W {v 1 - 1 )*< 



K*-i) 



/ 1 + B\* 
Observe that P is always associated with ( v 1 — ), and that, each form being 

repeated twice, there are two forms of function of each kind. Moreover, a cosine 
and a sine function are 'always associated with different kinds. It is obvious that the 
:ji -functions are expressible in terms of the ordinary P-functions of spherical 



harmonic analysis, and that if we take out the factor ( ~— — -=l£ ) the P-functions are 
similarly expressible. This factor will occur so frequently that 1 write 



a(p) 






l -/5 



\ V 



and as elsewhere commonly put O to denote ll(v). 
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We assume then the following forms for the functions : — 

For EEC, OEC, OOS, EOS 

$' = qf* + W% - 2, P ~ 2 " + 2/3% + .,„ P + *'. 
For EES, OES, OOC, EOC !-.... (6). 



P = 0{ ? ',P+ ^Y S .. 2 „P^-^ + 2£y. + 2.,P' 



,9 + 2-» ) 



In these series n proceeds by intervals of one at a time, beginning from a lower 
limit of unity. In both forms the upper limit of the first 2 is ^s or -J (s— 1) according 
as $ is even or odd ; and the upper limit of the second 2 is |r (i — 5) or -J (i — - s — 1) 
according as i and ,9 agree or do not agree in evenness or in oddness. 

The factor O contains (v 2 — • if in the denominator, but P* does not become 
infinite when v = ± 1, because when .9 is not zero P* is divisible by (i/ a — 1)* and 
we shall see that q is zero. # When s is zero there is no function of the P type. 

It may be noted that the limits of the series are such that neither q nor q' can ever 
have a negative suffix. 

We shall ultimately make q s and q s equal to unity, and this will be justifiable 
because there must be one arbitrary constant. 

We have now to consider the forms of the cosine and sine functions. They may 
be derived at once from the preceding results, for we have only to read (z/ 2 — 1)^ 

1 + 8 \^ 
as cos t<f) where t is even; (i> 3 —■ 1)* as sin<£, ( v 1 — * ~ •-- ~- - ] as cosc/>, and v as 

(1 — /3 cos 2<£)*. 

The factor (1 — /3cos 2<£)* will occur frequently, and 1 write 

®((f)) = (1 — /3 cos 2<£) ri , 

and as elsewhere I commonly write <X> to denote <3> (<f>). 
The following are the results : — 



% 



Type O or EEC ; <8? = $y t cos t<j> . 

.. 

v' ■ 

Type AB or EES ; 5* = S 7^ sin £<£ . 

2 

It is clear that we may equally well regard the lower limit in the latter as zero. 
Type A, or OOC ; each term is of type cos (t — 1) <j> cos <J> or cos (t — 2) <j> + cos £ <f>. 

i 

Hence (£,* = 2 y/ cos £<£ . 

1 

* This also follows from the fact that the series for P s begins with ttcc 2 (v 2 - 1) in the case of EES, 
and with £la s v (v 2 - 1) in the case of OES. Thus in the former case there is no term Oa and in the 
latter no term Qa\v 
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Type B ? or OOS : since we now have cos (t — i) <j> sin cf>. 



jj* = % y i sin /*</) . 



Type C, or GEO : 



O* = C P S 7/ cos (^ — .1 ) (jy 



' Type ABC, or OES ; each term is of type <$> cos (i — ])</>sm <f> cos <f>, which gives 
sin (t + 1) c/> — sin (t — 3) </>] <J>. Hence 



S 9 = <I> 2 yt sin (£ — I.) ^ . 

It is clear that we may equally well regard the lower limit as unity. 
Type CA, or EOC ; each term is of type <J> cos (t — 2) & cos <j>. Hence 

O = ®kyt cos (t — ].)<£. 
Type OB, or EOS ; each term is of type <E> con (< — 2) <f)sm<f). Hence 

S* = *Sy/sm(« — 1)0. 



When i and s agree as to evenness or oddness we have the forms independent of <[>, 



when they differ in this respect the factor <E> occurs. 

Therefore (in alternative form) for EEC, EES, OOC, OOS 



c 



= «.(«> +Sfe 



?'* 



2«< 



COS , 

• (•' 

sin N 



2w) </» + %fi"p i + 2n \ C0H (s + 2».) 



and for OEC. OES, EOC, EOS 



• (7). 



s 



= <n 



/ r cos ,,_*„, r cos , 

^ sin 



r co 



sin 



2^<£ + W. + ^ sil > + 2«)tf. 



In these series n proceeds by intervals of one at a time, beginning with unity. In 
both forms the upper limit of the first 2 is ^s or \ (s — 1) according as s is even or odd. 
In the first form the upper limit of the second 2 is | -,-(i — s), and in the second form 
it is \ (i — s — 1). 

We shall ultimately put _p 5 and ^/ y , which may be regarded as arbitrary constants, 
equal to unity. 



§ 5. Preparation for determination of the Functions. 

In order to determine the coefficients q, q', p, p> and cr, we have to substitute these 
assumed forms in the differential equations. 

Where the functions involve O and <I> as factors, the forms already given for the 
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differential equations are perhaps the most convenient, but in the other cases a 
reduction seems desirable. 

By considering the forms of 1) L and I) 3 in (3) it is easy to show that 



xli 



(**-l) 



L 

dv 



i(i + 1)(^ - 1) - * 






£ 



(^ _ i)( v 2 + i) -f- o+ 2^-f - t(t +1)^-0- 



rfr 



• • 



X* 



d~ 
defy* 



+ s 3 - /3 



COS 2* --,■; 



rZ 



sin 2<£ ~-~ -j- {(i + l) cos 2</> -j- cr 



> » 



(8), 



(9)- 



By making /3 vanish we reduce these operators to the forms appropriate to 
spheroidal harmonic analysis. By making /3 infinite we obtain the differential 
equations specified in § 2 as appropriate to ellipsoids of the class c 3 = -|-(cr + //). 

It is now necessary to perform the operation xjj 6 . on typical terms I 3/ and fl P*, and 

I COS I COS 

v y on typical terms < . ' t<b and <3>« . * £<A. 

/V J1 [sin ^ [sin ' 



(a.) To find i/r,(P). 

The form (8) for i// 9 is here convenient. 

It is clear that 



(" s - 1) 



dv 



i(i + 1) (v a - 1) - s'-j F = (z 2 - s 2 ) F, 



because P' is the solution, of the differential equation found by erasing the ter 
— s 2 P' from each side. 
Again we have from the same differential equation 



$ 



(" 9 -l)^P 



dv 2 



2v~ + i{i + 1)F + -j^--: P. 



rfp 



z/- 



1 



It may be noted in passing that this is equally true when the subject of operation 
is Qf 9 the function of the other form. 
Therefore 



dv 



» d 



(^-l)(^+l)£ + ^-.-(t + l)^ 



(Yz> 



cr 



I 



:>/ 



— 2i> 



+ ii + i) + « 2 - - 



<^£> 



Hence 



^(F) = (i 2 - s 2 )F - £ 



O 



^ 






-j- z ( ?/ -j- 1 ) -4- & 



1 



o V 2 -f 1 



cr 



_j 



F. 



F 



1 



O" 



i 
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e have now to eliminate i> — and ; 7 r. 

c/p v — 1 

It i. knowa that P = .,V,('f )V - 1)', 

and P^ = (/ - lfYf-YP a 

The differential equation satisfied by P^ involves £ in the form /A Hence 
(v % — l)"W-y- ) P can only differ from P^ by a constant factor. In order to find 
that factor suppose v to be infinitely large ; 

then 



and 

Also ( 

Therefore the factor is . — . and 

% — t\ 

I (1 V i 4- t\ < d\~ l 

F = („* _ ivWii p = l±-i^ („* - i)-*'(-f ) P. 

It will be convenient to pause here and obtain the corresponding formula} for the 
Q-functions. Various writers have adopted various conventions as to the factors 
involved in these functions. I write 





- ' >-"(£) 


P 

F = 

-i 


2i\ . 
2i\ v l 




9, 


2H\ i - tV 

9 7 ! -? t 


2i! j/* 




^(i!) 8 i + *! 


- 2 { i ! i + * i 



.CO 



vo' *— — ■ JL 



A*?- i)(P) 2 



and Q< = (* 3 - iW|;Yq. 

As in the case of P' we may change the sign of t 9 if we introduce a constant 
factor, and this may be found by making v infinitely great. In that case it is easy to 
show that 

2* (i if i 



Q 



2i 4- 1 ! v i + l 



d y _ / d \ 



By performing ( y- 1 and ( -j- j on Q it follows that the constant factor is the same 

as before, and that the alternative forms for Q* are exactly the same as for P f \ 

Hence the transformations which follow for the P-functions are equally applicable 
to the Q-functions. 
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If we differentiate P / in its two forms we find 



iW 



dv 



fdV /dV + 1 



dv, 



And 






i - *! 



£j> (i> 2 



dv 



1) - 1 „.„ ( lj' P+( , 



_j? p i —i 

y 2 - 1 ^ (v* - I)* 



1) 






dv 






tv 



v 



2 



i + f- t\ ' i~+ t - 1 1 <y - rj* ' 



I now write 



{i, t} = (i + t) (i- t + i) = i{i +i)-t{t- i). 



It is clear that 



{*, 



«} = {t, e + 1}, and {i, 0} = {i 1} = t>' + 1). 



i -\- t\ i t 4- 1 ! 

Now since -7 — 7-! — = \i, t\, by taking; the sum and difference of the two 

% — t ! ^ + t — 1 ! l J ,y ° 



forms of -y-, we have 



<£P 



dv 



vV f 



1 

2(>~ - 1) 

1 



(^ - 1)* 



L 

2 






-r 



£ 



p* + i _|_ 



p-i 



> 



(10). 



It is easy to verify, by means of the relationship P~~ l 



i~zJJ 



V\ that these 



equations are true when t is negative. They are also true when t = 0, although the 
second equation then becomes nugatory. 



v 



Multiply the first of (10) by v and the second by _ -1— ^ and apply them a second 



time. 

Then since 



K * + 1) 

^ + 1 

{U 4- 1} 



t - 1 



2i(i + 1) 
~ / 3 - 1 



4. '--* _ o/^i-.l)^ , 

*(* + 1) ""*"" *(* - 1) "~ ,\ *» _ 1 



I/ 3 + 1 



1 









Z>" — I 



i; 



2v 



clV 



dv 



v" 4- 1 

^ Z p* 

2 1 

i>" — 1 



9 



1 

2 



* + 1 

* p^+3 

/(TTi) 



* 3 -- 1 



+ ^ZT-"-'- F "* 



2^ (4 + l) p/ , {i 0{ij - 1} p/ _ 2 " 

f~ - 1 "•" *(1 - 1) 



> ■ • ( n )- 



These equations are always true although for i = ± 1 and they become 
nugatory. 



3 P 2 
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Then 



(I pt v 2 4- 1 
2^— + z 3 -~ 1^ 

«z/ zF — 1 



I 



F + 2 + {i, t] {% t - l}F-8] - i(i -f 1)F. 



Hence 



,/,, (F) = - i£ 



r\i ' p/ _j_ p/ + ? 



2<rF + {M}{M - 1}F 



(1.2) 



(A) To/me/ ^ (ft F). 

It is now best to use rj/, 9 in the form (5), where D t is defined by (2). 



Now Fj (OF) = ■ 



(i - m 



J v i — l)(l — yg) — 2>S 



a 



i\n 



i — ~t- — p^ 

civ v° — 1 



and 1)/ (QF) = (^ - 1) a [(z/ 3 - 1) (1 - /3) 

*-p ( 1 "^ p 

\ 
The latter terms of \p s contribute 



&J kJ 



,pj>i 



civ 1 



y , 



(3 \ 9 rfF 



2/3F /i/ s + 1 



z^ 



1 



v 



1 



12 { — i (£ -f 1) [(?/ — I) (1 ■— /3) — /3 1 P / — (,t — /Scr) P'} 



Therefore 
^ (X1F) 



cU 



(*f "~ 1) [(^ ~*~ 1) (1 "~ fi) 



2 



dir 



+ T(V 2 — 1) (1 — iS) + i8] 2v 



dV' 



civ 



i(i+ 1)(1 _- £)(„«_ ])p 



s 8 F 



~y" /j£ (^ 



But (^ 



1 ) - 



2v -y + {(/ + ;i)F + ■„ 

civ v- 



/ 3 



+ 1) F + 0<rP' - 2j8P' ^-J ► . 



P'. and we find on reduction that 



xji, (nF) == np- (* 3 - **) - ySO 






}/ + '.,-"" " (t 2 + 2) P ,; — cr 



F 



vcl L zF "4** I 

On substituting for ---: and -,,- .'■ P / their values, we have 

av v u •— 1 



B 



t/»,(ftP0 



:j8n 



r 2 F 



/ 2 ) 







F 



I 7^ " p/ + 2 O ~ p/ 

i 






t 



i > 



'j '^ '// j -j £ ? t *~ " I j x '' 



* s 



. (13). 



/ r ( \ 

( r .) Tofirulx.^y^ x t;t>J. 

In this case the most convenient form for ^ is that in (9), and we easily find 
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/' [cos , 
/V \ sin T 



P 



2(s 3 -f) 

y ~ 



GOS , . i • i -a i 

Sill 7 ; 



COS , , ^\ . 

• (* + 2)<A 

sin x ' T 



cos 



+ 2<r\ '." «A + {*,*} \\ (t 

•sm r ' c ) sin v 



1 



2)<f> V . . (1.4). 



1 now use the form ^ as defined in (5), where D a is given in (2), so that 



I) = (l — flcos 2 <i>)^ ', = <D 4v 



We have 



D Q (<H C0S *<A 

\ 1 sm '' 



3 



and 



sm f 



t « < si1 : up ± i/s (« + i) { ;:;: <« + 2) * ± # <« - 1> \ zu - 2 > * 



cos 



COR 



<fc 



9.\ J 008 



« s {^^ + ii8(f + l)(«+2)|^(«+2)^ 



+ ^(*-l)(«-2){^(*-2)^ 



The latter terms of ^ contribute 



* 



cos 



(*-Po)\*rti>-tfi(i + i) 



< cos (* + 2) </, + 



cos , * , 

. (t — 2)6 

sm x 7 ' 



rn 



Therefore 



/ *• cos , , 



PJ> 



2(*' 



/3 



" :) "G^ + &* + 2 HZ(' + 2 )4> 



cos 



sin 



+ 2cr\ . t6 + {i, f 



1 !{:;;:<'- 2 > 4 • (»»■ 



§ 6. Determination of the Coefficients in the Functions. 

In. this section I use successively the four results (12) (13) (14) (15) obtained in 
the last section under the headings (a). (/3), (y), (8), 



( a ) 



19* = 2 ,p + Zpq.^Y*- 2 * + S0^ + 2 /P 



+ 2/7. 



The limits of the first 2 are 1 to |s or h (s — 
or | (i — ,s — 1). 

Applying the operation \p s to SjjJ* and equating 



1), and of the second 1 to r> (^ — s) 



2 



~ ^ (p*) to zero, we have 
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is + 2n 



$8n(s - n)p~ l q,_J?-*» - %8n{s + n)^~ l q, + i J? 1 

+ #,[P +2 - 2o-F + {i, s] {i, s - 1}P~ 2 ] 

+ tpq, _ 2» [P ~~ 2n + 2 — 2crP s ~ z " + {i, s - 2w} K s — 2n 

+ $P"q, + m [P + 2H + 2 - 2crP ,+ - to -f {?; 5 + 2m} {{, ,s + 2n 



i } p - 2n - 2 ; 



0. 



The coefficients of the P ? s must vanish separately. This gives from the coefficients 
of V s " 2 " and P* + 2tt the following : — - 

2[4n(s — n) — /3ar~] g 9 „ 2n + /%^- 2 

+ {i, 5 — 2n + 2} {?', .9 — 2n + 1}^ _ 2 » + 2 = 0, 

— 2 [4n (s + w) + /3cr] g s + 2n + ^ + 2n „ 2 

+ /3 2 {i 3 ^ + 2n + 2} {i, 5 + 2n + l}q 9 + 2n + 2 = °- 

These equations may be written in the form 

2q s -2n ~{i s ~ 2n -\- 2} {i, s — 2n + 1} 

*7* — 2n + 2 j / \ n i i n9 /^$« — 2n — 2\ 5 

4% (5 — w) — po- + f p^ ( — - 1 



9 



<?s + 2 



% 



#s - 2u / 
1 



$s + 2% — 2 



4n(s + w) + /So- — i/3 3 {?', s + 27?, + 2} {^, s -f 2% + 1} ( 



/2ft 



.? + 2w + 2 



2* + 2 



to 



. (16). 



Whence by continued application, the continued fractions 



2q,. 



s — 2n 



{% s — 2n 4- 2} {i, s — 2?? ~f 1} 



: /3Hi s - 2n} {i, s - 2n ~ 1} 



\ 



<Z - 2n + 2 



4?i(s — n) — f3cr — 



4 (n 4- 1) (s 



w 



i) 



i/3 3 {is - 2» - 2'/- + 2}{i,3- 2w - 2r + 1} 

ff* - 2n - 2r - 2\ J 
Qs - 2ti - 2r 



4 (?i + r) (s — ^ — r) -f ^/S 3 ( 



9, 



^.s + 2% 
f /.s + 2n ~ 2 



1 



4?i (s + w) 4- /So- — 



i/8 2 {i, s + 2^ -f 2} {i, s + 2n + 1} 

4(^ + 1) (s + ^ + 1) + /So- — ... 



(16). 



J/3 2 {i s + &+ 2r} {i s + 2^jf 2r --J^ 



4(ti + r)(s -f % + r) - i)8 2 {i,s 4- 2n + 2r + 2}{i«s + 2n + 2r -f i}/^±^± 2 ^ +2 \ 



We must now consider what I may call the middle of the series, which corre- 
sponds with n = 0. In this case each of the t's contributes one term, and the g s term 
mves another. The result is 

- 2crq, + fiq t - 2 + Pq 9 + 2 {% s + 2} {% $ + 1} = 0, 



or 



/3or = 1/3* 



9, 



9.8 -2 



+ i/3 2 K ^ + 2} [i, s + 1} 



'2j 



l± 2 



Since 2q s ^ 2 /q s and 2g, + 2 Ai, are expressible as continued fractions, we have an 
equation for /3cr, if the continued fractions terminate. 
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We shall now consider those terminations. 

First, suppose that s is even, corresponding to types EEC, GEO. 
The first continued fraction depends only on the first 2. The condition to be 
satisfied is 

+ 0*? O [P» - 2crP + {i, 0} {a, - 1}P-*] 

+ j8* ( - 2) g 2 [P* - 2crP 3 + {% 2}{i, 1}P] 

+ j8^- 4 >^[P 6 - 2crP± + {-/, 4}{i, 3}P 2 ] + ... = 0. 

Since {i, 0} {i, — 1} P ~ 2 = P 2 , we have, by equating to zero the coefficients of 
P and P 2 , results which may be written 

2 ?o _ -fi2tKl | 2ft _ - Ji4}{j3} 

Hence the g's disappear from the first continued fraction, which terminates with 

- iiS 8 {i 2} R 1 } 
s 2 — ^jo- 
in this last term the ^fi* which prevails elsewhere is replaced by |-/3 3 . 
Observe that when 5 — 2 the first continued fraction is replaced by a simple 
fraction, so that the equation for /3cr becomes 

/3a = Hi^i^ill + ^ R 4} {i, 3} (&) . 

Secondly, suppose that 6' is odd, corresponding to the types OOS, EOS. 
The condition to be satisfied is now 

2 ( S 2 _ l)£*<-»> ?1 pi + 2(s a- - 9)£'<- s >2 3 P 3 + ... 

+ ^-'^[P 8 - 2CTP 1 + ft 1} {i, Oj P- 1 ] 

+ i8* ( -^ 8 [P 5 - 2o-P 3 + ft 3} ft 2} P l ] 

+ P*-<»< h \F - 2crF + ft 5} ft 4} P»] + ... = 

Now ft 1 } ft 0} P " '= i (i + 1) I 11: P - i = i (i + 1) P 1 , and if we equate to 

ii J. ', 

zero the coefficients of P 1 and P 3 we obtain results which may be written 

2fi __ nl*'lL^^l 

<? 3 ' S 3 - i _ y9«r + i/3* (i + 1)' 

■i'h rJA. 5 'f,l? ; ' 4 ] 

?6 s 3 - 9 - /3,r + iW 2i M 

V 'i3 / 



o. 
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Thus the </b again disappear, and the fiist continued fraction ends with 



,3 _ i _ fa + X^ + !)• 



Observe that when s = 3, the continued fraction, reduces to a simple fraction, and 
the equation for /3<x becomes 



P<r 



— 1R* {') 9 >i li *>\ l : hi " 

i£ J±>.:L>J± ".; . . l. i^i^ (jh; 4ip s 



The case of s = I must be considered separately. 

At/ 

We have next to consider the termination of the second fraction, which depends 
onlv on the second S, 

First, when i and s are either both even or both odd, the types are EEC and OOS, 
and the limits are |- (i — 5) to 1. The condition to be satisfied is 



+ 0*('-*- 4 >#_ 4 [F- 2 - 2dP ; - 4 + [ii - 4}{i,i - 5}P ; -°] + ... 



0. 



Now P' + s is zero, and equating the coefficients of P ; and P' 2 to zero we obtain 
results which mav be written 






I 



tf,„ 3 ^ - s a 4- £cr J 






(i _ 2)2 - ^ + £0- 



1 

|/3 3 {i^}{i^l}(^~) 



Hence this continued fraction, ends with 



I- 3 — s 1 4 /3cr 

Secondly, when i and s differ as to evenness or oddness, tlie types are OEG and 
EOS, and the limits are ^ (i— s— 1) to 1. The same investigation applies again when 
i is changed into i — - I. 

Hence the continued fraction ends with 

-.jfi*{^- 1}{m- 2} ^ 

(^ — l) 2 — s 3 4- @o' 



The cases of ,9 = 0, s = 1 must be considered by themselves. 
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When s = 0, the types are EEC and OEC. The " middle " of the series is now 
also an end, and the condition is 

- 8&F - 8 . 2 3 y % t P i - . . . + g [P 3 - 2crP + {i 0} {i, - 1}P~ 2 ] 

+ Al 3 [P J - 2o-P 3 + {i, 2}{i 1}P] + ^[P 6 - 2<rP 4 + {», 4}{i, 3}P e ] + ... = 0. 



Writing P 3 for {i, 0} {i, 
zero, we have 

f$ar 



1 } P a and equating the coefficients of P and P 3 to 



Si 



^{i,l}{i,2}Q 






1 



4 + per 



mi 3} k 4} (^ 



Therefore /Jo- 



1/3 3 K1}K2} 
4.1 s + /So- — 



i/3 2 {i3}R4[ 
4 . 2 3 -f /3<r - 



—.1/32 k a f^ ^ _ -| I 
ending with —*—^-r^^ — for EEC, and with 



i* -f y8<J 



4.3 3 + /3<r - ... ' 

} F{i,i- l}{i ,i- 2 
(i - l) 9 + /9<r 



for OEC 



Next when s = 1 the types are OOS, EOS ; the " middle " is again an end, and 
the condition is 



-8.1. 2g 3 P 3 -8.2. Zpq.J?* - . . . + q x [P 3 - 2oF + {i, 1} {i, OjP" 1 ] 

+ /8<7 8 [P S - 2crP 3 + {t, 3} {*', 2JP 1 ] + /3% [F - 2<rP 3 + ft 5} {1, 4}P 3 ] + ... = 0. 

Writing i (i -f 1) P 1 for {1, 1} {i, 0} P ~ l and equating to zero the coefficients 
of P 1 and P 3 , we have 

fa. - y3i(i + 1) = -l£ 3 K 3 }{», 2} f^ , 



\2i/ 



2ft 

ft 



Therefore 

y8cr- ££4(1 + 1) 



i9o--i/8»{*,5}K4}(^) 



4.1.2 + 



jffi 3} {j, 2} 
4 . 1 . 2 + /3<r - 



jff te 5} ft 4} 
4.2.3 ~f /3cr - 



i^K 7} {*, 6} 



4.3.4 + ftcr — 



ending with 
EOS. 



?' 3 + /8cr 



for OOS, and with 



£/9 3 {i,i- 1}K* - 21 



(* — 1)- + /3ar 



for 



(A) We have next to consider the other form of P-function for types EES, OES. 
OO C, EGG, namely, 



p^ 



where 



,2 _ 1 ±J\\ 



a 



J,, 



1 - /3 



IP 



1 
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Let us write q' t ± 2n = (s i 2n) q t ± 2rt . The </s are not now the actual coefficients 
of any P-function, hut we shall see that they are determinable by almost the same 
relationships as those already found, and therefore the notation is convenient. 

We now have 



6 



n[g,*P' + tp%„ ail (s - 2n)P"" 2 " + tp% + 2n (s + 2w)P + aw : 



Applying the operation i//? to P 5 and equating 



<> 



0/3 



t/^. ( P*) to zero, we have 



$8n(s - n)(s - 2n) q,_ 2n P*- 2n t8n (s + n) (s + 2n) q 8 + J? 8 + 2it 



+ 9.s[(s + 2)P + 2 - 2crsP* + \i 9 s} {i, $ - 1} (s - 2)P 



s -2" 



+ 2j8* 3 , _*[(*- 2* + 2)P-» 

+ {% s • 



n. + 2 



s - 2a 



2n\ u, «s* 



2cr(s — 2w) 

— 2?i — • I j Is 



2n 



2 



\ p« - 2». 



+ 2j8«g, + 2n |> + 2w + 2) F 



+ 2h + 2 



2cr (5 + 2n) P + 2n 



+ {i, ,9 + 2n} {i 9 s + 2n — 1} (s + 2n — 2) P 



a 4- 2n - 2*" 



0. 



This is the same equation as before, if we replace tP l by 'PK As we may equate 
coefficients of iP* to zero (instead of coefficients of P*), we obtain, the same equations 
for the q's as before. 

A certain change must, however, be noted with respect to the beginning of the 
first series, which determines the end. of the first continued fraction. 

We previously wrote V 2 for {i 9 0} {i 9 — 1} P"^ 2 and i (i + 1) P 1 for {i 9 1} {i 9 OjP"" 1 . 
But the corresponding terms will now be {% 0} {i, — 1} (— 2)P"~ 2 and \i 9 1} {i 9 0} 
(— 1) P" 1 , and these are equal to — (2P 2 ) and — (1 . P l ) s 

Hence it follows that when s is even (EES, OES) 






{i, 2} {% 1 

3 2 — /3<T 



<?4 



S 3 



-- 4 — /9<r " 



The </ term has disappeared from the latter of these, and thus the continued 
Taction is independent of q . This is correct, since whatever value (short of infinity) 
q may have q' 09 being equal to 0q , vanishes, Hence the continued fraction is 
docked of one term and ends with 







4. — ' £jQ~ 



It is important to note the deficiency of one term in the fraction,. since it indicates 
that when s = 2 the first continued fraction entirely disappears. 

When s = there is no function; of the P form, so the question of interpretation 
does not arise. 
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4 8 M 



When s is odd (OOO, EOC) the only change is that i(i +1) enters with the 
opposite sign, so that the first fraction ends with 



l&{i, 3} {% 2} 



s< 



1 — /3cr 



- %&i(i + 1)* 



When s = 1, we have ficr + g/&(i +1) equal to the same fraction as before. 

When the g's are determined we have q t = ^. But it is desired that in the 
case (a) q s should be unity, and that in the case (j3) q f s should be unity. This con- 
dition will be satisfied in the present case if we determine all the r/s, put q s equal to 
unity, and finally take 

s 4- 2n 



C l s ± 2n 



( l$ ± 2n 



Thus in both (a) and (/J) we put q 8 equal to unity, and in (/3) determine the q's by 
the above equation. 



(y.) We now have to consider the cosine and sine functions. 
For EEC, EES, OOG, OOS 



©* 
& 



sm r L 



cos . 

• ($ 

sm x 



r pos 



The first £ has limits |-s or Jr (s — - 1) to 1, the second -J (i — 5) to 1. 

2 / f OTA 
Apply the operation x* & n d equate — ^ x* (1 a* ) to zero > then 



J bill [ bill 



+ jp# 



{{, 5 + 1} -I , * (s + 2) <b + 2<x 4 . $<4 + (1 «s} 1 . (s — 2) 6 

1 > 1 j 1 gm v 1 /r 1 ^ gm r 1 1. j J [sm v 7 \ 



T~ ^P P - 2n 



cos 



cos 



{i, g _ 2n +1} . (s - 2n + 2) <£ + 2cr 4 (5 - 2w) <£ 



sm 



sm 



[cos 



+ {i,s - 2n} j sin (5 - 2w - 2)<£ 



+ ZP"p 9 + 2n 



{i> s -f- 2n + 1} * 



cos 



. (5 4- 2n + 2) d> + 2cr - 

sm v - / r 



fcos 



. (s + 2n^ <& 

sm x ' ^ 



+ {i, s+ 2n\ \ T (s + 2n — 2) <£ 
1 J [sm x 

3 Q 2 



= 0. 
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cos 



If we equate to zero the coefficients of i J - K (s 4- 2n) <b, we find 



s — 2» 



{% s — 2n + 2} 



P*-2»+2 4 7l ( fS _ ft,) _ ^ _ JL£2 {^ s _ 2 



W 



-7^ +_2« 



-~{i } s + 27?/ — 1} 



1!- (~P s -- n ~ 2 



('% 



4n ( s + ,,) + fie + 1/33 {,; * + 27i + 2} ( ^*± 2|, + ? ) * 



These will, as before, lead to continued fractions, and by elimination of the ps to 
an equation for /3cr. The equation will agree with our former result, for it can of 
course make no difference from which equation we determine cr, # It follows then by 
comparison with the previous result (16) that 

1 



Ps — 2w 

Ps + in 
Ps + 2» - 2 



{i, s — 2w + 1} &_ 2 „ + s 



{&, 5 -f* 2?i — 1} 



2* + 2» 
?•* -f 2k - 2 



Hence when the (/s are found, the p's follow at once, 



(8.) For OEC, OES ? EOC, EOS 



c 
s 



cb 



p/ < . St 

1 ' sin 



cos 



r<£ + 2/3^ _ 2}l <! T (5 — 2n) (f) + X/^/, + 2h. 

j bill 

where <S> = (1 — /3 cos 2<£) 1 . 

The limits of the first S are Js or ^ (s — 1) to 1, of the second \ (i 
Proceeding exactly as before we find 

2p f s~~9n {h S ~ 2ll + 1} 

__ - - -• • - 

' *P *-2»-2 \ 

. _____ 1 



. (s + 2n) <p 

Sill v 7 7 



_»*-s«+2 4^( s _ ft) — ^Qj- — |-/9 3 {^ s — 2w} 



{i, s + 2?4 



^P * + 2» ____ _ ___ 

Z' + sn-a 4^ ( s + ft) + ysv + J/3 3 {i s + 2w + 1} 



^P ,y _f.g 



2%-b 



i 9 * + 2« 



By comparison with (16) we see that 



P *-2» 

JP «?- 2m + 2 



1 



5* -2m 

K s — 2;i + 2} ^^ 2 „. F3 



— {£, 5 + 2n} 



PVf2«~2 ' _ / * + 2»-2 

Therefore when the ^ ? s are found, the p s follow at once. 



1) to 1. 



* I have of course verified that this Is bo. 
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We may now summarise our results, as follows :— 

In the general case where s is neither nor 1, /3or is the root which nearly vanishes 
of the equation 

o ~ i& ft 4 ft * - 1} i$* ft * - 2 } ft s - 3 } 

pa =■ 



+ 



4 . 1 (s - 1) - /3(x - 4 . 2 (s ~ 2) - /3<r — .. 

i/3 3 {i s + i} ^ s + 2} i/3 3 ft s + 3} ft s + 4} 



4 . 1 (s + 1) + /3a- - 4.2 (s + 2) 4- /3<x 



The continued fractions terminate variously for the various types of function, The 
end of the first continued fraction is as follows : — 

— i/9 2 ft 1} ft 2\ 

For EEC — ~ — — — -—- ; and when s = 2 this is the whole fraction. 

5" — per 

— i/3 3 ft 3} ft, 4} 

For EES — * — -~ — - L — ; and when 5 = 2 the fraction disappears. 

For OOC T~^7^~^hr^ — 7x I an d when s = 3 this is the whole fraction. 

s - — 1 — . fie — •fpft'?, 4- 1) 

_ i/93 17 2} ft ?A 

For OOS ^~~^~~-^^~-^^-^ T; ; and when s = 3 this is the whole fraction. 

s -i __ J — ftp 4- -|£fft7, 4. 1) 

— x /9 2 {i 1} ft 2} 

For OEC —^-5 — '—z — " — J and when s = 2 this is the whole fraction. 

S~ — per 

— -i-iS 3 U 3} K 4} 
For OES ■~~*^~ - ^r~^ I ~~" ; an $ when 5 = 2 the fraction disappears. 






— 1 /3 3 ft 21 It 3} 
For EOC ^-^r~ J ~7r~ L:L T7^T' — T7 5 an( l when 5 = 3 this is the whole fraction. 

r — 1 —• per — f/&(t + 1) 



I/53 ft 21 ft 3} 

For EOS - v~ ± ~i^~~H L T^7'~T~^ ? an d when 5 = 3 this is the whole fraction. 

S" — 1 — per 4- fpi (^ 4 1) 

For the first four of these types, viz., EEC, EES, OOC, OOS, the second continued 
fraction ends with 

— ~:r""~o"7~£~~ ~ J an ^ when 5 = ?! this is the whole fraction, but with the 

%" — S" 4" per ' 

sign changed. 
For the last four, viz., OEC, OES, EOC, EOS, it ends with 

X S 2 ft, i — - 1} ft, i — - 2} 

71 — 77, — ~r~i^w~~-- 1 and when s = i — 1 this is the whole fraction, but 

(7, — 1)- — &"* 4- pa ? 

with the sign changed. 
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When ,9 = 0, the equation becomes 

ifi* {i, 1} \i, 2} i^K3}{i4} 



fia 



ending when i is even (EEC) with 



4.1 3 + fia - 4.2 3 + £<r - ... ? 



i° 4- /3cr 
- i/3 2 K i - 1} K i - 2} 



and when i is odd (OEC) with ,. iV> 

x 7 (t — 1.)'- -f /3cr 

When 5=1 the equation has two forms, which may, however, be written together. 
If the upper sign refers to cosines (OOC, EOC) and the lower to sines (OOS, EOS), 
the equations are :— 

i/&{i,2}{i3} t/3* {i, Q {i, 5} 



ficr ± \fii{i + 1) = ;. 



ending when i is even (EOC, EOS) with 



4.1.2 + ficr -— 4.2,3 + j3(r 



(i - If - 1 + fi<r 



when i odd (OOU, OOb) with --■•— v~r~5 • 

It might appear at first sight that a difficulty will arise in the interpretation of 
these results when i is small, for the numbers in the denominators of the fractions 
increase, and yet it is possible that the number at the end should be smaller than that 
at the beginning; thus apparently the fraction ends before it begins. But this 
difficulty does not really arise, because in such cases the numerator will always be 
found to vanish, and thus the whole fraction disappears. For example, in the last case 
specified, if s = 1, i = 2 the denominators, according to the formula, begin with 
8 + ficr and end with + ficr ; but the fraction has for numerator {2, 2} {2, 3}, 
which vanishes. 

When ficr has been determined we find the q J s by the formulae - 

2q.~2n - {i s ~~ 2n + 2} \i, s - 2n -f- 1} \^{i y s - 2n} {i, s - ^-JL} 



q t -2n + 2 ^ n ( s — n ) — (So- — 4:(n + 1)0 — n ~~ 1) — £<r - 

2q, + 2n 1 \j&{% ■&_ + 2ti + 2} {i,? +_2n__+ 1} 

( h + 2n~-2 ~~ 4^0* + ™) + /3cr — 4(ft + l)(s 4- ^ + 1) + fi<? — ... 

The terminations of the continued fractions are as specified above in the equation 

for ficr. 

By forming continued products of ratios of successive </s, we can find all the q's as 

multiples of q s> and q$ = 1. 

In the cases EEC, OEC, OOS, EOS, these are the required coefficients for $*. 
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c» —i— 2i')"i 

In the cases EES, OES, OOO, EOS we put q! ±2n = - z= ; — ^ ± 2rt5 and thus find the 

coefficients for P*. 

The coefficients for C, & in EEC, EES, OOC, COS are determined by 

px ^_2rt_ _ 1_ ^* :l 2h_ 

i^« - 2» + 2 {'^ s 2?l -f 1} q s -2n +2 

P* + 2n i • rt , -i l ( ]s 4- 2?i 



{i, $ + 2rc — 1} 



P* + 2w - 2 3* 4- 2 



n 



The coefficients for G, S in OEC, OES, EOC, EOS are determined by 

Ps - 2n _ _ j- . ,ll~ 2 !L 

p'* - 2n + 2 ~~ {is- 2?1 + 2} q a __ 2n + 2 ' 

_£« + **_ _ _ ^ 5 + 2n } -A±^-... 

P s + 2w - 2 ' "?« + 2/i - 2 

It follows that if we put q s = 1 and _p, = 1 



_ , , i 

P»~2n — ( — f {^TT2^Tl}K« ^Ito, + 3}T. . {i,*~- 1}* 



s - 2n ) 



jft + 3 » = ( — )*& 9 + 2n — 1} {i, s.+ 2n, — 3} . . . {i, s + 1} sr, + 3a , 

jp'« + sn = ( — )* {'' : > * + 2 ?i} {i s + 2w — 2} . . . {i, s + 2} q, + 2v . 

When s = 0, Qg/go * s equal to that which would be given by the general formula for 
Jlif±i ? i when we put in it n = 1, s = 0. Hence it follows that the </s for 5=0 have 

9.9 +2» — 2 

double the values given by the general formula. 



If we change the sign of s, the two continued fractions in the equation for ficr are 

simply interchanged. Hence /3cr is unchanged when s changes sign. Also, since 

{i, t} is equal to {•— i — ■ 1, t], ficr is unchanged when — i — 1 is written for i. A 

consideration of the forms of the q's and p's shows that q^ s + 2 k ~P~ 8 + 2k is equal to 

'' — si 

. — -7<7„-2fc P**" 2 *, and therefore 

7,- -~p »s ; 









P/~ IP'-i-, 



488 PROFESSOR G. H. DARWIN ON ELLIPSOIDAL HARMONIC ANALYSIS, 

■ § 7. Rigorous determination of the Functions of the second degree. 

If a numerical value be attributed to (3 it is obviously possible to obtain the 
rigorous expressions for the several functions. Thus, if /? were J- we could determine 
the harmonics of the ellipsoids of the class c 2 = \ (a 2 + 6 2 ). But I do not think it is 
possible to obtain rigorous results in algebraic form when i is greater than 3. In 
order, however, to show how our formulae lead to the required result I will determine 
the five functions corresponding to i = 2, but I will not work out the case of i = 3, 
although it is easy to do so. 

When s = 

— M!] 2 ^ ^ I?i 2 ! — 12 ff * 

4 + /3cr 4 -f /3<r* 

Therefore /3cr == -—2 + 2(1 + 3/3 2 ) 1 , or writing i> 2 = 1 + 3/3 2 for brevity, 
/3cr = 2(B — I). Then putting gb = 1, and remembering that the value of q% is 
twice that given by the general formula, 

1 B - 1 

J2 — 4 + /Jo- - " ~6£» 

Therefore |3 3 = P 3 -| — — P 2 2 ... (17), 

w 



here P 2 = |^ 2 - £, P/ = 3 (^ - 1). 

The coefficient of: the cosine function is given by 

2h = — ( 2 > :L ) ( h = — G 2:^ 
Therefore C 3 = 1 — — tj— cos2</j . . . . „ , . . . (18). 

5=1, cosines ; EOC type. 

The continued fraction has ^f$ 2 {i> 2}{{, 3} in the numerator, and vanishes because 
[2, 3} = 0. Therefore 



/3a- + }j/3i (i + 1) = 0, where i = 2 . 
Therefore or = — 3 . 

But the coefficient is independent of /3, for 

Ps 1 = O fei'JV]* and ?/ = 1 . 

Therefore P, 1 - A/!L_i^.p o i == 3i/ i; 3 - f^£ . (19). 

Clearly CV = cos <f) (1 — /3 cos 2</>f . , (20). 
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s = 1, sine ; EOS type. 

The continued fraction again vanishes and a .= 3, but it is not needed to express 
the functions. Putting q( = 1, 

W = q [ 'V l ^? 1 '=3v(S-lf (21), 

S 2 ] = sin ^(1 - ft cos 2<f>)* . (22). 

s = 2, cosine ; EEC type. 

The second continued fraction vanishes because it contains {2, 3} in the numerator. 
The equation is then 

_ -$/F{2,l}j2,'2} _ -J2/3 3 
P(r ~~ I - £<r " 4 - '/So- ' 

Therefore /3<r = 2 - 2 (l + 3j8 a )* == 2(1 - J?). Then putting q, — 1, 

„ _ _ i {2i 2 H_ 2 Li^ _ - K^jz.i) 

io ~ ^ 4 - /So ~ /3 3 
Therefore f a 2 = — - ( | -— >- P 3 + P 3 3 (23), 



where P. 



an -I 



2 



o^ 



2? 



P/ = 3 (i/ 3 - 1). 

1 
Then p = - -^-y q = - £«/ , and 

&/ = ^- + cos 2<£ (24). 

5=2, sine ; EES type. 

Both fractions disappear and <x vanishes, but is not needed for determining the 
functions. Noting that q ' = 0, and q% = 1, 

P s s = n[ % 'P^ = 3^-^|y^-i)* (25), 

5S () "zzz Slu Zi(p * . . . . . « . . , , . , , , , ( ZiK} ). 



OV' 



We can write down the functions of /x by symmetry, and the products of the three 
functions give rigorously the five solid harmonic solutions of Laplace's equation of 
the second degree. As I have remarked above, the seven harmonics of the third 
degree may be obtained rigorously by a parallel process. 

§ 8. Approximate Form of the Functions. 

It is clear that the first approximation to /3cr is zero, and that the second approxi- 
mation, in the general case, is - . . 
VOL. cxcvii. — a. 3 R 
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y8cr 



I rz ii $} {% f — 1} , i pa {i> s + !} fc 5 + 2} 



s 



1 



P 3 (35 s - 2i(i + 1) 



s + 1 

^ a (^ 4- 1) 2V 



5- 



If this expression were inserted in — - — we should obtain q s± % correct to /3 s . But 

since the next approximation would only introduce /3\ it follows that q s ± % would be 
correct to /3 3 inclusive. Now q s±2 enters in the functions with a factor /?, and there- 
fore this approximation would give results correct to /3 h inclusive. Since the similar 
operation could be applied with equal ease in all the cases in which the continued 
fractions assume special forms, it follows that this degree of accuracy is very easily 
attainable. However, the forms of the coefficients would be rather complicated, and 
it would render the subsequent algebra so tedious that I do not propose at present to 
carry the approximation beyond /3 2 . 

It now suffices to put a = in the denominators of all the continued fractions, 
whereby the coefficients are determined, except in the cases of s = 1, s = 3, where 
we put cr = ± ^i(i + 1). 

In the general case we have 



( h-i 



1 {i s} {j sj ~ 1} 



8 



9* -4 — 



2*-2 



tf*-* 



{i, s}_ {i, g — 1} {i , s — 2} \i, s 
128(5 - r)(g^2) ~ 



s 



9, 



s — 4 



2*-2> 



2*--4j 



P*-l 



\i 8} 



'}>*-* 



1 J s-% 



P S-4: 



8 (s ~ 1) \ 
{% s} {i, s - 2} 



0\ 5 



128(5 - 1) (s - 2) 
{*', s — 1} 

{js — 1} ft, s -~ 3} 
~128(s - l)(s - 2) ' 



3} 



* + :2 



f i^ + 4 



f i<? + 2 



5 *+ 4 



P^-f-3 



P* + 4 



P * + 2 



P^ + 4 



1 

8 (s 4- 1) 



\ 



1 



128 (s + 1) (s + 2) 

s 4- 2 



s 2*-f-2>' 




s 4- 4 

O - -■ - 




— {i s + 1} 




.8(s 4~ 1) ? 




{i, s + 1} K 6- 4- 


3} 


128 (s-4- l)(s 4- 


2) ? 


- {i,5 + 2} 




8(5 4- 1) ? 




{is 4 2} {i, s 4- 


4} 



128 (s + l)(.s + 2) J 



■ W- 



When s = 0, we double the results given by the general formula and find 



% 



There are no 



4 > ( li — "li?> .Pa 

t/,', ^,', and p. z ' 






. . (28), 
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When s = 1, 

q. A — _ ,, — rr = Te(l ± it A' (i +1)), with upper sign for cosines 



(EOC, OOC) and lower sign for sines (DOS, EOS). 
9h =: T2"8Ti!T3 =: TtF 8" * or a *f cases. 



But for OOS, EOS we use the P form, and for "EOC, OOC the P form; and for 

the latter = 3, = 5. 

s " s 

Therefore for OOS, EOS (sines) 



and for EOC, OOC (cosines) ... . (29). 

<7s' = A (1 + -iV # ^ + 1) ), q 6 ' = y|- 8 - ! 

For OOC, OOS, with upper sign for cosine and lower sign for sine, 

Ps=-i i s{i,2}[l±Tk&(i+l)], Ps = vhs {i 2} & *} ■ -(29), 
For EOC, EOS, with upper sign for cosine and lower sign for sine, 

Pz = - i\- K 3 \[l ± -A-fr (i + I ) ] , p-: = jh K 3} K 5 } . . (29): 

When s = 2 the coefficients may be derived from the general formula. 
When 5=3 

fc = mh^H) = ~ A *'"' 2 > {i 3 ' f [1 ± lW * + 1)] ' 

the upper sign applying to cosines (OOC, EOC) the lower to sines (OOS, EOS) ; 

9.5 = 32 5 9l == 2T60' 

But for OOS, EOS the |S form applies, and for OOC, EOC the P form applies. 

s — 2 ___ x s -f- 2 __ 5 § 4- 4 



Also with 5=3, = i 



s - 3 ' s — 3 > s ~ * 



3 R 2 



492 



PROFESSOR G. PL DARWIN ON ELLIPSOIDAL HARMONIC ANALYSIS. 



Therefore for OOS, EOS 



Qi — ~ -h fr 2 l {*'. 3} [i - iV^' + i)], g 5 = A. ? 7 = 25V0 • v 



For 00C, EOC 



( h 



-is{h2}{i,3}[i+-^^(i+l)],q : / 



9 6 ? r i7 



__7 
7680 • 



For 000, OOS, with upper sign for cosine and lower for sine, 



Pi 



iV {*', 3} [1 ± -&{ti (i + 1)] , p B = - ^ -{/, 4} , 



i^7 = ¥5lToK4}{^6} . 



For EOC, EOS, with upper sign for cosine and lower for sine, 



P/ = iVK2}[1±^/3^+ i)LiV 



5 2 I ^ * J i ? 



1 } 7 — !J560l i J J itM) 



(30). 



It will save much trouble to note that if we were to admit negative suffixes to the 
q's, the general formula would give us the term /3^q^ l P"™ 1 , where 



2-1 



{i3}{i,2}{i t l}{i,0} 
128 ."2.1 



1 



Thus this term is —r tt /3H(i + 1) . [i, 3} {i, 2}P l 6 But this is exactly that part of 

the term in (30) which arises from fiqj? 1 , but which is not included in the general 
formula. 

Similarly the general formula gives for q f ^ v P-i> P-i those parts of the terms 
arising from q/, 2 } \> Pi which are not included in the general formula. 

It follows that in much .of the subsequent work we need not devote special 
consideration to the case of s = 3, 



§ 9, Factors of Transformation between the i too forms of P function and 

C~ or Sf unction. 

The rigorous expressions \$ s and P* always differ from one another, but approxi- 
mately they are the same up to a certain power of /3 ? provided that s is greater than 
a certain quantity. 

2/3 



Since fl = 



2 __ l+j\k 
1-3 



= 1 + 



2 . 7T3TT77~q\ ) > it is legitimate to develop fit in powers 

of l/(p 2 — 1) up to a certain power, say t, provided that it is to be multiplied by a 
function involving at least (V 2 — l) f as a factor; for this condition insures that there 
shall be no infinite terms when v = Jr 1. At present, I limit the development to /3 3 , 

so that 

j3 -4- /3 3 li& 






1 
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Therefore 



— - -*&*—) p* + / 1 







\ 



v- 



v (i8 9 ',_ 2 p- a + ^'. + a P + 8 ) 
+ /6 Y. - 4 P " 4 + A 



9 ' 

/ 8 + 4 



ps + 4 



It is obvious on inspection that we cannot rely on this development if s is less 
than 4. 

If then s is equal to, or greater than 4, this value of P*, when properly developed, 
to the adopted order of approximation can only differ from H^ s by a constant factor, 
say C* or shortly C s ; so that 



f * = O P 



(31), 



and we have to determine the constant G\ 

We might develop the above expression for P* completely and compare it with 
P*, but this is unnecessary since the comparison of a single term suffices. 

I now write 

*._ *(* + !) 

2,1 ~ ' s 2 - 1 ' 



(32), 



or shortly % This notation is introduced because this function occurs very frequently 
hereafter. 

We have seen in (11) (slightly modified) that 



P» 



P» 



-f o 



. i . , _ 2 (S + l) P + I^i^r ii p-2 ^ . 

zr — 1 [s(s 4- 1)- . 7 s(s — 1) 



We may write this 



I 



^s 



V" 



a,P + 2 + /8,F + 1 + y,p- 2 , 



where 



a, =- 



Then 



45 (s + 1) ' 



P' 



ft = — H^ + ] )> r* 



{?; 4 {^g — i} 

"~ 4s (s - 1) 



(i/ 3 - I) 3 



= a,(a, + 2 P + * + & + 2 P + 2 + y. + 2 P) 

+ A(a,P + 2 + AP + yj P- 2 ) 

+ y,(a_ 2 P + /3 s _ 2 p- 2 + y_ 2 p- 4 ). 
Therefore the coefficient of P s is a 8 y s + 2 + (/3,) 3 + a, _.. 2 y„ or 



16 



s(« + l) 3 (s + 2) "•" ^ "*" ' "*" s(s - 1)* (s - 2)_ 
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I now introduce a further abridgement and write 



J. t 



(i _ i) j (j + i) a + 2 ) 

s 3 - 4 



• • » 



or shortly T. 

Then, after reduction, I find 



(32), 



P* 

<y - Y) 3 



1[- 2V + 2 s + 42 + 3 + TIP + . . 



Accordingly the coefficient of P* in P* is 

1 + ^(2 + l) + i/3 3 (2 + 1) - tV^C- 2 2 * 3 + 2 3 + 42 + 3 4- T 



But g\< 



i/3 3 


</s - 2 




L(« - l ) ( s - 


-2) 


(s — 2) {-j, «} {% * 


— 1 



+ q's + 2 



{a, s + 1} {i, g + 2}' 
(s + l)(s~+ 2)~ 



$ 



q.* 



s + 2 



? 7 s + 2 



and the last term in the 



8s (s - 1) ? ** + * 8s (s + 1) 

above expression will be found to be equal to + i@ 2 (t 2 ' — 1). Thus the coefficient 
of P* in the development of P* is 

1 + 1)8(2 + 1) + i\/3 3 (2V + £ 2 + 4S + 3 - T) ; 
but the same coefficient in 31* is unity. 

■p §y 

Therefore 



1 
C? 

l 

(C ; ? 



= 1 + &8 (2 + 1 ) + iV/S 2 (2 V + S 3 + 42 + 3 - T) 
= 1 - $/3 (2 + 1) + ■&/&(- 2V + 32 3 + 42 + 1 + T) 

= l + /3(2 + 1) + !/3 3 (2V + 32 3 + 82 + 5 - T) 
= 1 - /?(2 + 1) + |/3 3 (~ 2V + 52 3 + 82 + 3 + T) 



I . . (33). 



The squares of this constant and of its reciprocal are given because they will be needed 
at a later stage. 



We next consider the cosine and sine functions. 



c< 



= $ 



cos j. i a ' [cos, . , , 

. s<£ + ftp s _ 2 \ . (s — 2)6 + ... 

sin r ■* sm v ' T 



As far as f3r 



* = (1 - ficos26f = 1 - 4/3cos2<£ - -^(l 4- cos 4<£). 
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Therefore 



c 



1 _ %j3 cos 26 - te/3 1 (1 + cos4<^)] 



cos , 
. Sep 
sin 



4- R 






1 — |/3 cos 2(f] 



cos 



1 Sill 



C0,S / I A \ ± 

. Is + 4) 0. 

sm v } T 



This expression, when developed, must lead to (&? or 3t s multiplied by a constant 
factor. 



Let 



(JP* j\ s JO 

S 1 " M S 



(34). 



cos 



Then D/ or D* may be found by considering only the coefficient of < . s<f). Hence 

j bill 
1 



But 



V 



Therefore 



\y 



1 1 i$P 4H P 8 ~ 2 4,H V s + 2* 



{is 



11 



8 (s - 1) ' 



o 



p & 



4- 



• {^ s + 2} 
8(s + 1) 



P« - 2 + P'a + 



J 

4 



; (S + 2). 






i 



(D/) 3 



1 _^(S + 3) 



and 



• » » 



• . . « { o o ), 



The reciprocals may clearly be written down at once. 

There are no factors by which |3 3 , \$' 2 , p 1 can be converted into P 3 , P 2 , P l ; but 
this is not true of the cosine and sine functions. 

In the case of s = 3, it will be found that the general formula holds good for the 
factor whereby (£ 3 , j& 3 are convertible into C 3 , S 3 . 



When s = 2, 



C2 

S 3 



J;8cos 2c£ — iVi8 a (l + cos 4<£)] 



cos 
sin 



. n 2«£ + /3p' (1 ~ ££ cos 2e£) , 



+' /8p' 4 (1 - |/8 cos 2</>) . 4<£ + ftp 



lhen 



cos . . . ^ , cos 



Sill ' • ■ - D slu r 



C 3 = [1 - (-& + |p'o + ip\) /3 3 ] cos 2^> + . . . 
S- = [1 — ( sV + ii>'+) J ] s i n 2(b + ■ • ■ 



496 



PROFESSOR G. H. DARWIN ON ELLIPSOIDAL HARMONIC ANALYSIS. 



J3Ut p — -g- 1 1, 1 ; , J> i 



2 



nr {*", 4} , p'g = Tir Vo {i, 4} {i, 6}, and 



•32 + i^'i, = — st (2 - 5) , where 2 = * 9 3~-' ■ 



1 



|. ( '£, "T" L ) , 



Therefore the factors are 



1 

V 

1 
i) 



(cos)^ 1-^(52* + 7), 
i(sin) - l+A^(S 2 -5) 



(36). 



It is easy to verify that the other coefficients of C 3 and JS 2 are in fact reproduced. 

The notation adopted here and below for distinguishing the cosine and sine factors 
is perhaps rather clumsy, but I have not thought it worth while to take distinctive 
symbols for the factors in these cases, because they will not be of frequent occurrence. 
When s = 1, 



[C 1 






1 — |/3 cos 2cjf> — y^/3^ (1 + cos 4(j&) 



cos 



r* + #/ s (l-*j8oos2${^3* 



I /~jo / { COS _ » 



1T^-^(P'3 + 1)] {"> + 



This must be equal to T A ^ . 
Now, with upper sign for cosine and lower for sine, 



p 



3 



lVK3}[l ±-^i(i+ 1)], / 5 = T kK3}{^5}. 



Substituting for _p' s its values, we find with the upper sign 



1 



-(cos) = 1 



IP 



W (p' s + i) = i - i£ + -iiP \i{i + l) 



10 



"> 



And with the lower sign 

\ (sin) = 1 + i/S — i£ 3 (p' s + i) = 1 + P + Ti V/S 2 [t (i + 1) - 10] 



(37). 



I).. 1 



It follows that 



/ 



D/ (cos) 
I)/ (sin) 

"Di 1 (cos) 



i + &-■&? [»(* + 1) 



14" 



•2 



l-i^ + ^^p ({+!)- 8] 



= l+ii8-Jj/3»[t(i+l)-16] 



. 



* • 






l "l 3 

— (sin) = ! + ££ + ^ 2 [i (i + 1) 



D 



IV (sin)] 3 = 1 - p - -^[^i + 1) 



8] 
16 j / 



• ( 3 ?) 
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We cannot in the present case use 2/ as an abridged notation, because it is infinite 
as involving s* — 1 in the denominator. 

It is easy to verify that the other coefficients of (& l and % l are, in fact, reproduced 
in the transformation. 

Lastly when s = 0, we have only cosine functions. As before 

C = l — i/3 cos 2<£ — -£ 6 p 2 (1 + cos 4<£) + fa'* ( I —.1/3 cos 2<j>) cos 2f/> 

+ /3 2 9/ 4 cos £(f> 

= 1 - &/& - i£y s + • ■ • 

GC 1 

This must be equal to n , and therefore — - = 1 — y^/J 3 — ^i^p,,. 

Now jf>' 2 = - f[/, 2} , j>\ - T |s {i, 2} {i, 4} . 

Hence ■ = 1 + -j 1 ^ [»' + 1) — 3] (38). 



Since in this case %; 



i(i + 1.) 



- 1 
1 

I); 



— — I. """*"*■ •% 1\Kj £h i ) O. . • • ♦ • • a • i lOOli 



Thus the general formula again holds good. 

It is easy to verify that the other coefficients of (£ are in fact reproduced. 
The principal use of the transforming factors, determined in this section, is that it 
will enable us to avoid some tedious analysis hereafter. 



§ 10. The Functions of the Second Kind. 

The second continued fraction of § 6 terminates because 

{i 9 s + 2n + 2} {/, s + 2n+l}=0 
when 11= I? (i—s) or \ ({— •$— -1), since one of the two factors then assumes the form 

| '£, 2, —J- 1. j . 

Hence it follows that the equation for determining or is the same as before ; but we 
cannot on that account assume that the q coefficients vanish when their suffixes are 
greater than L 

In considering the P-functions it was immaterial whether or not we regarded them 
as vanishing, because P^ vanishes if t is greater than L But the Q-functions do not 
vanish in this case, and therefore we must postulate the existence of qs with suffix 
greatgr than i. 

VOL. CXCVII.— A. 8 S 
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In fact, whilst we have as before, when i and s are both odd or both even, 



2q i 



1 



( 1 



i~z 



% 



;z 



s 3 -f /5a- ' 



we also have 



2ft 



+ 2 



Qi 



(i + 2f - s» + /3a- - W{% i + 4} {% i + 3} 



' ^ft + 4 



V?- 



* + 2 



2#; 



and similarly a fraction for — - - 4 , and so forth. 

#' + 2 

It follows therefore that while the </s with suffixes less than or equal to i depend 
on finite continued fractions, those with suffixes greater than i depend on infinite 
continued fractions. 

It thus appears that while the first series in the expression for ©./ or for Q/ has 
limits 1 to j^s or J (s — 1), as before, the limits of the second series are 1 to qo . 

Thus we have found, an expansion for this class of functions in powers of /3. 

In the limited case in which the coefficients have been actually evaluated, namely, 
where the development is only carried as far as the squares of /3, we have 

1 1 

£* + * — 12S(g + i)(^~2) ' 



s + 2 



2 s+{2 



8s (s -f 1) 



, s + 4 

q s + 4 "" i28s"(s" + ' 1) "(s +T) ' 



These coefficients do not vanish when s + 2 or s + 4 are greater than i, and this 
confirms the conclusion already arrived at. 

In spherical harmonic analysis there is no occasion to consider the value of Q. s when 
s is greater than i 9 and the values are therefore not familiar. I will therefore now 
determine them. 

It is known # that 



Q 2 



2*(i \y 



2i + 1 



v 



1 

i + 



1 \ ( 



i + 2! 



1 



2 .1H! ' (2i + 3)* 1 



+15 + 7>o 



i + 4 ! 






2 3 .2!i! (24 + 3) (2i 4- 5) i> 



i + 



5 i* » * • 



Therefore differentiatin 



dv 



q. = (-y + iy.i! 



And 



I)'* 2 * 



t + 1 1 (i + l)(i + 2) 1 

<Ua»3a«M ~ " """ * "~~" i aaa t-nn _ .. .. . „ — , __ __ — .......... ._„ — ' 1 _ _ 



1 

.2* + 2 



I>" 



/ y + i 



2* . i ! 



<y. _ i)<+i' 



, , 2* + * . i 4- 1 ! 
(-) 1 — ^m:v. 



(v» - 1) 



i + 2 



* Bryan, < Camb, Phil 3oe, Proc.,' vol. 6> 1888, p* 203, 
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d V + 8 
dv 



Q I -=.(-)«2 i + 1 .*+l! 



2 (/i -f 2) z/ 






1 



_^ (;/ _ l)i+8 ' (^3 _ 1) 



\i + 2 



/ \i + 1 



2 l l + 2 .i + 2l 2 i + 1 (2i+ 3).i + 1 

\i + 3 » 



> 3 - iy 



(* 3 - I) 



* + 2 



/ d 



\ i + a 



\dvj 



But 



Qi = (-Y 



2 i + s ,i + 3l 2 i + z (2i+ 3).i + 2! 



> 3 - 1) 



i + 4 






(„2 - l)i + 3 



J>. 



Q< = (v s _ !)1* f^YQ., 



therefore 



0/ + 1 



/ y + i 



^ 



2'i! 



(„2 — 1)H* + i; 



q: 



i + 2 



Q 



i + 3 






> 



/ tooct \i + 



Qi* 



£ + 4 



(-Yh 






>i + 4 + (2i + 3) (i/ 3 — I)] 



(„3 - l)i(« + 4) 



[2i + 6 + (2i + 3)(v 3 - 1)] 



/ 



. (39). 



These are all the functions which can be needed for the expression as far as /3 3 of 
(&i or of Q/ when s is less or equal to ?'. If s is equal to i, we shall have terms 
ft* qi + 4 , Q s+4 or 0/3 2 g'^ + 4 Q* + 4 , and these are the furthest. 

But it is well known that there is another expression for these functions of the 
second kind. 

The differential equation is 



d 



(v 2 — l) 2 yi; + 2v(v" — 1) 



7 O 

dv 



d 
dv 



/3 



(ir — 1) (v + 1) 



(l/ 3 


— 


1) 





& 


ri 3 


+ 


2p s 


d 

dp 


. — 



^H-i 



i (i + 1 ) v% "~ or 



where (©,/' may be interpreted as meaning also Q/. 
Let us assume that 



Ydv 



m = o, 



is a solution, where <£§/, :}| s may be interpreted as meaning also Q*, P' 
Then since ffi is a solution of the differential equation, we have 



(v 2 - l) 3 



2V * + ^ ^ 



rf*/ 



eZi/ 



/3 



■+ 2i/ {f - l)p s V 



(^ _ i)(„» + l)(2V * + $»?) + 2,^'V 



<??/ 



<&> 



3 s 2 



0. 
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This is easily reducible to 



d 

dv 



log 



fc« 



v(f^)V - in* 2 



1 4" /3^ ''■ 



0. 



whence V 



Hence 



m i 



'if 



V 






F- 6/ 



where (3*/ Is a constant. 



i -p 



»x> 



©,* = e'^i 



r?z/ 



""•\ 



r(ftW-l)V-l^' 



Q/ = E/P* J 



» 



«.// 



F 



„(P/)V- i)V-{4f)* 



The general solution of the differential equation must be 



a$* + y 



i ? 



• • 



. (40). 



and we have already found both ffi and ©; ? . Hence the two <El* % & must be different 
expressions for the same thing, for the form of (@/ as a series negatives the hypothesis 
that it involves p* in the form yffi + y 3 ®*. 

Having then two forms of (Q? or of Q*, it remains to evaluate the coefficients (Bf> 
E/ 5 which are involved in the equations (40). In order to do this it will suffice to 
consider the case where v is very great, so that 



P^ 



2i\ 



2 i ir 



V" 

i - t 



q* = (-)%-; 



2* A\ i ± t\ 
2^+1! ~t^+ t ' 



i 



As far as concerns the first term in the series 



w = 



211 



v 



2H\ i - si 



1 + /% 



i — s I 



% — s + 



7; + #7* 



+ 2 



- s — 



9 



+ P\-a~; 



% — s 



— 5 + 4! 



; + /% 



+ 4 



^ 



i — s 



4! 



©' 






2\ ^ ! i + s ! 



?; + i! z^ +1 



i 4. * __ 2 ' ^ 4- s 4- 2 » 

1 JL Rn L.TJL — ' JL. Rn ' 



**«3 



7/ + S ! 



■s + 2 



'?.' -j~ 5 I 



+ A*?-4 i+jJ ; +^^ + 4--+^ 



It will be observed that if s is equal to 1 or i — 1 the terms in ?j§f in g, + 2 and gw 4 
disappear; and if 5 is equal to i — 2 or i — .3 that in gw 4 disappears. This agrees, 
as it should do, with the vanishing of V s + 3 and P* + 4 when the order is greater than 



the degree. 
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7 



It we write ?$ s = av { and (©* = y+v ^ ie ^ rs ^ °f our equations (40) becomes, when 



v is very large, 



7 



z; 



* + i 



= (Q* ay* 



tCa d 



v 



, v v* (otv l f 



<B 6 



V' 



.00 



d 



V 



& i v v 



Q.j* 4- Q 



31 + 



®' 1 

a v 






i + l 2% 4- 1 



Therefore C* = (2i + 1) ay, and since the a, y in the case of the P*, Q s only differ 
from these in the accenting of the qs we have 



e= (-Y-. 



i + s I 



% — 5 1 



i + ^_ 87 



i — s! 



i — 5 4- 2 ! 



I * * • 



11/3 1 -{• S — L \ 

1 -]- pq s _ 2 7TT^ + • • • 



-i 4- s I 



E* 9 = the same with accented q's. 
Effecting the multiplication of the series 



<B? = 



x i 4- si 
(-)' — 



£ — 51 



1+/3 ?,_■ 



^ — s ! 



i - s + 2 ! 



"T Qtf-2 



-i + s — • 2 



+ g 



^ 



5 4- 2 



% 



s\ , i 4- s + 2 V 

- « i... . /-y 



5-21 



i + 5! 



+ P(q. 



% — 5 I 



5-4 



i — s + 4! 



+ ?*- 



i -f s — ■ 4 ! 
i -\- si 



+ 5', 



5 4-4 



^ — 5 I 

i — 5 — 4 ! 






5 + 4 



^ -f* 5 4~ 4 I 

i + s! 



+ q*-*q*-* 



5 ! i -f 5 — 2 ! 

5 ! i — 5 4- 2 ! 



' i — s\ i + s -h 2 ! j — 5! i + 5 4- 2 ! 



j- 



2, 4- 5 I ^ — s — ^ 1 



5 ! i 4- 5 — 2 f 



E* = the same with accented g's. 

If we substitute for the qs their values, the coefficient of /3 inside [ 



in the 



expression for (ffiE* is 



J, 

8 



(i 4- s)(i 4- s — 1) (i — 5 -4- l)(i 



5—1 



5 — 1 



5 4- 2 ) (^-^(i-g-l) 






5 4- 1 

(?' + s + l)(i + s + 2)' 



5 4- 1 



In the expression for E s the first pair of these terms are multiplied by 

5 4-2 



5 - 2 



and 



s 



the second pair by 



s 
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The coefficient of /3 3 in the expression for (& s is 



128 



+ - v -- -'- v -- -r 1 —^—— — + 



(i + s)(i + s— l)(i + s- 2)(i +j — 3) (i -_s + l)(i-s + 2) (i - s + ?>) (i - s + 4) 

_...___ j_ ^ ^ ^ ^ _ 

(i-«)(»-8-l)(*-s-2)(i-s-3) (i + s + 'i)(i + 8 + 2)(i + 8 + 3)(i + s + 4) 

-r ,. . u,, , ox ~r (s+ i)( s + 2 ) 

2(i - s )(i - s - l)(i + s + l)(i + g + 2) 

(s + I) 3 

2(i - sj- l)(i-s)(i-s+ l)(i- s + 2)1 
~~ (a - 1) (,s + 1) ~~ ~~ (8 - 1) (sTl) j " 

s — 4 
In the expression for E* the first pair of these terms are multiplied by ; the 

s + 4 , . As — 2 x3 i i A + 2-\ 2 
second pair by ; the first of the third pair by ( ) , and the second by f 

a* - 4 









(* 


+ i)( 


s + 


2) 








2(i 


+ 


S )(i 


-{- ,9 - 


-i)(* 


— 5 


+ 


i)(i 


— 


8 + 2) 








(s — - 


i) 3 










2(* 


+ 


s — 


i)(i 


+ s)(i 


+ s 


+ 


i)(i 


+ 


s + 2) 



and the last pair by 

o 

The reduction of terms such as these will occur frequently hereafter, and I will 
therefore say a word on the most convenient way of carrying it out. It is obvious 
that the coefficient of (3 may be arranged in the form 

M (i + 1) + B (2i + 1) + C . 

The coefficient A is equal to the coefficient of $ in the original expression, and if 
we put i = 0we have B + C, and with i = — 1, — B + C. Hence A, B, C may be 
easily determined. 

Again the coefficient of /3 2 may be arranged in the form 

A# (i + I) 2 + B (2i + 1) i (i + 1) + Gi (i + 1) + D (2i + 1) + E . 

This may be written 

Ai 4 + 2 (A + B)i 3 + (A + 3B + G)i 2 + (B + C + D)i + D + E. 

It is easy to pick out the coefficients of £ 4 , r\ i 2 , and we thus obtain A, B, C. 
Then putting i successively equal to and — 1 we have D + E and — ■ D + E. 
In order to express the results succinctly I use as before the notation 

V - - i l± 1 ) T . _ (* ~ l)*ff+l )(t+2) 



s 



2 1' * s 3 -~4 



and I usually omit the superscript and subscript $ and i. 



PROFESSOR G. H. DARWIN ON ELLIPSOIDAL HAEMONIC ANALYSIS. 503 

Proceeding in this way I find 

m = (~Y r^, (1 - MP - 1) + A/8»[- *"(S 9 + 22-1) 



« — s i 



+ 3 (2 2 - 2$ + 2) + 2T] } 

(-^TTTi ^ + W& + 3 ) + liW* 8 ^* -22+1) 

- (^ ~ 26S - 42) - 2T] } 



. (41). 



1 
These results may be verified, for if we multiply (& s by --—, as given in (33), we ought 

to find E* ; and this is so. 

The formulae apparently fail when 6* ='. 0, 1, 2, 3; but when s = 3 they still 
hold good because, as remarked above, the general formula for s = 3 gives correct 
results when properly interpreted. Thus it only remains to consider s = 0, 1, 2. 



When 5=2 the coefficients of /3 remain as in (41). In the coefficients of /3 2 

9s_4 = 0? 2* -2 == *§" I 7 "? ^j (.&> 1; ? 9* + 2 " 8"T3 ? 2* + 4 " X"2 8"Tl3"T¥ 

?*--4 == ^, ?*-2 == , tys + 2 = TT5 •> ?* + 4 " T 2*8"71T • 

In the expression for ^H 2 the coefficient of /3 2 inside the bracket is 

nsk? [3 (i - 2) (i - 3) (i - 4) (i - 5) + 3 (i + 3) (i + 4) (i + 5) (i + 6) 

+ 72 (i — 1) i (i + 1) (i + 2) + 8 (i — 2) (t - 3) (^ + 3) (t + 4) 
- 24 (i + 1) (i + 2) (i + 3) (i + 4) - 24 (i - 3) (i - 2) (* - 1) i] . (42). 

Effecting the reduction and writing S for j$i (i +1)? w ^ find 

<&* = \±\\ { 1 _ p (2 - 1) + ¥ |^ 2 (19S 2 - 1302 + 80)} . . (43). 

The coefficient of /3 s for E 2 may be got from (42) thus : — Multiply the first and 
second terms by 3, erase the third, fifth, and sixth terms, and multiply the fourth 
term by 4. 

Effecting the reduction we find 

E! = i^j\ {1 + J>S(S + 3) + Th/8*(25S 2 + 1862 + 368)} . . (44). 

Observe that there is no factor by which p 2 can be converted into P 2 , so that this 
case cannot be verified like the general one. 
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When .s' = 1 we have 



Qs - i 



0, 



<7« 



0, 



q\ _ 4 = 0, (/, -2 = 0, 5 



(7* -i- 2 



_! 

16 L 



(7« + ^ 



•s* + 



3 
"16 



1 --K^(i + 1)] ? 



768 
76 8 



The terms in /Jgw g and /3q' 8 + 2 now contribute to the terms in /3 2 . 
For C5 l the term in /3 inside the bracket is 

¥.*[(* " 1)(* - 2) + (i + 2)(i + 3)] = {-[i(i + 1) + 4]. 



The term in /}*, of which the first portion is carried over from the term in /?, is 






- 2) + (i + 2) (i 4- 3)" 

2) (i - 3) (i — 4) + (i + 2) ({ + 3) (i + 4) (i + 5) 

4. 3(4 - i)(i _ 2 )(i.+ 2)(i+ 3)]. 



This is equal to — T -|-g- [r(i + l) a — 56i(i + I) — 180]. 

As we cannot now use the abridged notation with %\ which is infinite, I write 






ius 



e 1 



i + 1! 



* 



1! 



j = i [i -|- I). 
-_1 + p(y + 4 ) - yk/3 e (i 3 ~ 56/ - 180); 



. (45). 



For E^ 1 the coefficient of /J is three times as great as before, and the coefficient of 
ff l is 

T fs?<>' + 1)0' (i + 1) + 4] + tts 1 ^ [ 5 (* — !)(* ^ 2 )(* — 3)(i - 4) 

+ 5 (i + 2) (i + 3) (i + 4) (i + 5) + 27 (z - 1) (i - 2) (i + 2) (i + 3)\ 

On effecting; the reduction I find 



Immhu .7 ' 



-t + 1 i 



> + f /3 (./ + 4) + jhF^bf + 376/ + 1044)]. . (46). 



When s = we have only (& to determine. Here 



2.-4 = £ 



0. 



(/• 



« + 



J 

4? 



The term in ft is ^[i(i — 1 ) + (i + l) (i + 2 ) 
That in /3 :3 is 



<1* + 4 



iO' + i)- 



__1 _ 
12 8* 



lk[>'(* ™ 1)(» ~ 2)(i - 3) + (i + l)(i + 2)(i + 3)(* + 4) 

+ 8 i (i _ i) (i + i) (j + 2 )] - -A- (5/ + 14;' + 12). 
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Fherefore 



n^ 



<&i = i + ijQO" + i) + jkP(5j* + uj + 12) 

= 1 - 1/3(2; - 1) + «5W52 a _ us + 12) 
since 2» = — i (i + 1) = — j. 

Collecting results from (41),. (43), (44), (45), 46), and (47), 

(* > 2) m = (-yl^ju - ^(s _ i) 

+ TTf/^C- « a (^ 3 + 2S - 1)+ 3 (2 s ~ 2S + 2) 4- 2T 
(ffif? = ~\ {1 - £0(2 - 1) + ^6/8 2 [192) 3 - 130S+ 80] K 

' «R l = ~ H~T! ^ + *^0" + 4 ) - TkW - % - 180] K 
<& = 1 + i/50" + 1) + AjSW + 14/ + 12). 

+ 3V/8 3 1> 3 (3S 2 - 2S + 1) - (2 3 - 262 - 42) - 2T]f, 
E? = l , ± 4 y \ \l + 1/3(2 + 3) + 2k/3 3 [252 2 + 1862 + 368]}, 

E; = _ l±i.j {1 + ^y + 4) + ^5^(55^ + 3 (. 6> - + 1044)j] 

where t = •?-+--» T = <! -ili^tjid+i), y = , •(,-+!) 



. (47). 



(48). 



PART II. 

Application of Ellipsoidal Harmonic Analysis. 
§ 11. 2%e Potential of an harmonic deformation of an Ellipsoid. 

A solid harmonic, or solution of Laplace's equation, is the product of two 
P-functions of v and of fi respectively, and of a cosine or sine function of <j>. A 
surface harmonic is a P-function of /x multiplied by a cosine or sine function of <f>. 
'0 We found 

$'(„) = F(„) + tj3%_ 2n p-*»(v) + 2j8*£ i + JP- + 8 »(v), 

where P* (p) = — — rp f — - J (*/ 2 — 1)' ; and a similar formula held for P*(v). 

Hitherto we have supposed P' (/x) to have exactly the same form as Y l (v). But 
since /x is less than unity this introduces an imaginary factor when t is odd, and 
vol. cxcvii. — a. 3 T 
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makes the succession of P's alternately positive and negative when t is even. As 
this is practically inconvenient I now define 

p w = S# (£) <* s - :| >"- 

and then retaining the former meaning for the q coefficients, we give the following 
definition—- 

with a similar formula for P* (/x) e 

Thus we need only remark that in the functions of fx the qs corresponding to odd 
powers of /3 enter with the opposite sign from that which holds in the functions of v, 
and the whole of our preceding results are true with this definition of P^ (/x). 



If v defines the ellipsoid to which the surface harmonic applies, we require the 
expression for the perpendicular p on the tangent plane at v , /x ? </>, and that for an 
element of area of the surface of the ellipsoid at the same point. 

By the usual formula 



O O 



7 O 

*" *" i_ V_ , 

f v ("o 9 - ^J) 3 * 3 (v - i) 3 "*" kW 

(1 - 0) (jfi - l±&) ^_l . ,^(1-/3 cos 2<f>) 

(tf - m») (tf - x 7-T*) (49) 



Let dn, dm, df be the three elements of the orthogonal arcs corresponding to 
variations of v, /x, cj> respectively. 

Then by the formula at the end of § 1, 



...._ ? 



^ ^j ivw - 1) (,/ - ;-^|) / 

V Wdp) *» (1 - ^) (K±J - ^) 

df V ^o 2 - L ^i i ^ i ) (^T^r* - ^ 3 ) 



. (50). 



M^|* ^ - e™f (1-/3 cos 2<£) 



Therefore ( ~ = ^ LrJ^A !.-?. ^/i 1 ...,.. P J . . . (50), 

U'#/ 1 - /3 cos 2<£ v ' 
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and 

/ pdmdf \ 



v o{ v * ~ l)f^ a 



o 



d p d 

dn Jc 2 V(flv Q 



1 + ff 

1 - A 



(1-/3) 



/ 1 - B cos 2$ o\ : 



\ 



(1-/3 cos 2<£) (1 - /* 3 ) (i±| - ^ 2 ) 



I (50). 



Two functions, written in alternative form, 



"^^"Ipr Z 



V 



L*V 



m 



v< 







V 



»/(/*) <&M , 



are solutions of Laplace's equation, and together form a function V continuous at 
the surface of the ellipsoid v = i/ . Reading the upper line we have a function 
always finite inside the ellipsoid, and reading the lower line one always finite 
outside. Hence V is the potential of a layer of surface density on the ellipsoid v , 



and by Poisson's equation that density is equal to 



4t7T 



dY dV 

(outside) — — (mside) 



dn 



dn 



Let the surface density, which it is our object to find, be 

a surface harmonic multiplied by the perpendicular on to the tangent plane and by 
a quantity p. 



T-. » Cv *p ci> 

an hr v dv 



j\. 



P 



4:7rJl 



$V r 



%' ("o) i; *U ("o) - ©*'("<>) £■ W to) 



r^, 



o 



<rfe 



o 



But 



€/ (v ) = &/ m (v ) 



oc 



rfy 



^0 



m (»)t (^ - 1)* {v* - i^) j • 



Differentiating this logarithmically we find 



4:7tJc 2 P q (Vtf — 1) S (^ 3 



x _j- «\i j a constant}. 



1 - j3/ 



Hence surface density ^p/(^)C/(<j&) • P> where p is constant, gives rise to potential 



f inside 



4.7r/i 



{outside theelli P S0id ffi; v-o 



"^(tf-iJW 



z; 



££Mw*' 



i! 







j^ 



f / (/*) <£/ (<£) 

» * • i It/lJi 

The same investigation holds good with $/(</>), or with P, Q, C, S in place of 
the corresponding letters above. 

O 1. /j 
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Imagine that the surface of a homogeneous ellipsoid of density p, defined by v , 
receives a normal displacement Sn, such that 

Sn = p . ef / (jl) C/ (<f>) . 

Then the equivalent surface density is_p. ep'ffli (j^)C/(^>) 5 and we can at once write 
down the expressions for the internal and external potentials by means of (51). 

If x Q , y , z Q be the co-ordinates of a point on the surface, it is clear that the 
co-ordinates of the corresponding point on the deformed surface are 



/ A$/ ( M ) ©/ (<£)\ / /~6f/(/,)©/(^. 



/y „- /v | 1 | 



U \ 



o \ - » ^3 



Hence the equation to the deformed surface is 



9. . o '2 



9 9 9 

or since P7^" l i + 1) + &i (/_ i~) "*" /fcV ~ 1 ' ^ may be written 

(^ - "o 3 ) 5 = 2ef ,' 0*) C/ (<f>) . 
If we substitute for ~ its value from (49), this may be written in the form 

Jr 

This is the equation in elliptic co-ordinates to the deformed surface, but in actual 
computation the form involving rectangular co-ordinates might perhaps be more 
convenient. 



§ 12. The Potential of a homogeneous solid Ellipsoid, 

It is well known that the potential of a solid ellipsoid externally is equal to that 
of a " focaloid " shell of the same mass coincident with its external surface. 
If p f be the density of the shell defined by v Q and v + Sv, we have 



iir¥p f 



v* + 2^ - fr|)W + 2 ^ - l)*(*o* + 2^)* 



v* - — ±— Y (v^ - I)^ = inkPp ( V - t^q! W - l)^o - 



1-/3, 



V° i - fiJ 
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Therefore 






or 



p f v 8v = 



OV 1-/3+1-/3 



• ♦ 



• « » 



(53). 



If Sw be the thickness of the shell at the point where p is the perpendicular on 
the tangent plane, 

So ii! 
n = ^aO^ 



0' 



.P 



y»2 

If we multiply both sides of (53) by ~, we see that the surface density of the 



focaloid shell is 



V 



P^W" !)(V -f=-!) ^ 



OZ/ n l - js ^ 1 - ^ 



|8/ /v 

. — * — 







F 



«>. • 



k 2 



If therefore we can express _ 

x pi 



in the form of surface harmonics, it will be easy 

to write down the external potential of the ellipsoid by means of the formula (51). 
Before doing this I will, however, take one other step. 
It is easy to see that 



3v 



o 



V , L+J 8 _ of 2 _ 2 + s 

1 ~£ + 1 - f3 V ° ''3 (1 - £), 



^ 



2 







2 — i? 
3<T~/3)y 



. . . (54), 



where for brevity B = (1 + 3/3 3 )*. 

Now on referring to § 7, (17) and (23), we see that 

¥» W = P 8 (v) + ^j~ ■- P 9 a (v) , » 9 » (v) - - 



2ji? -_1) 



P 2 (v) + Vi (v) 



where 



p 2 (") 



- 3, .2 



2 



V 



1 

2 ? 



If then we put 



w (*) 






it is clear that 



a 



i? - 1 + S^S 



a 



and 



2 

a 



2/3 


? 


O ( jO ~~~ X ~" 


-£) 







J5-2 





3 (1 - /3) 



P 8 ! W = 3 (v« - 1) 






r 



r 



7 



— J3 ~f* 1 — /3 . 

20 ' 

._ b + 1 + 30 . 

-5- 2 



«' 3(1-/3) 
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3 
It is obvious then that our expression (54) is equal to — - -^ft (v ) ;p 2 2 (v ). Then 

3 1 -f B 
since — — 7 = — -- , we have the surface density of the focaloid given by 

w KprP — — _— Jf£— — I -- "0 -- - «L 3_„ 






But since z^T~1t+5\ + ^rzrr\ + z^i = * ( 55 )> 

With the object of writing this function in surface harmonics, and besides to 
enable us to express a rotation potential in similar form, we have to reduce x 2 > ?/, z 2 
in the required manner. 

I now drop the suffix zero, since we are not concerned with any particular 
ellipsoid. 

Referring again to § 7, (18) and (24), we have 

&* W = 1 - ~^ cos 2<f>, <&*{<$>) = —p + cos 2<f>. 

If then we put 

2(1- I?) „ b - 1 + p 

' — 9 /' — 77 — 1 — 3/3 ^ 

op 

we may write 

C 2 (<£) = ecos 2 ^ + 4 <&*(<f>) = e 'cos 2 <£ + £'• 

Let us assume, if possible, 

^^=4±|) = ^*» 0*) & 2 W + W M «V (</>) + #, 

or (> - J-^) cos* ^ = F(ap* + r )(ecos 3 «£ + £) 

+ (r(oLfi 2 + y)(e'cos 2 (f) + £') + ^ 
From which it follows that 

Fat, + GaX = 0, i? y £ + OyX + JT = 0, 

jpae + (jae 7 = 1, Fye + CryV = — - 

1 — p 
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rhese equations give 

F = 



1 



1 



e e' 



G 



i c 



a » ( ~ ? 



H=- 



a 

e 



a' 



'f 3 



i c 



and the condition 



ye 



7 e 



(i-0(3-$+<i + vt . 



0. 



N 



ow 



76 



1 -f J3 - 5 

~ r~ "is ' 

- 1 - 3/3 + B 



\ "T+/3 



76 

r 



1 + /3 + i? 
1-/9 ' 

- 1 - 3/3 - 



i? 



1 + /3 



Since these values satisfy the condition amongst the coefficients, the assumed form 
for aP is justifiable. 
I find then 



F- 

Whence 



1 + B B - 2/3 



2B • 3(1 - py 



G= - 



1 + B B + : 



AB 3 (1 - /3)' 



if = 



1 + /3 
3 (1 - /3/ 



3« 2 



#(»•-£!) 



1 + * * ^ »a(/0<M*) + ^ 4^ »sW «.*(*) + 1 • (57) 



2£ 1 + /3 



4B 1 + /3 



This is the required expression for x 2 in surface harmonics. 



Next assume 



- V 



m^&M + GiW m ®z 2 ($) + Rv 



If we put 
we have 



& 3 3 - 1) 

<&, (<f>) = e x sin 2 <\> + t lt <&£ (4) = e/ sin 3 <f> + & ; 



2f/? - 



Q g- i) 

2, 



*i = 



«i'=- 



- g + 1 -f /3 

/3 
i? - 1 + 3/3 



and 



3/3 



(^ - 1) sin 3 <£ = i>^ 2 + y) (e x sin 3 <£+&) + ^i( a > S + O (*,' sin 2 <£ + £/) + H v 

Whence F l9 G l9 H l have the same forms as before, and the condition to be satisfied 
by the coefficients is 

2.- 1 - — y^± -l & __ -SL — o 
*Si a Si Si Si 

It will be found that the condition is satisfied, and that 



* T i = 



l + s 



65 



) J 



^ __ 1 +J 



jfiTi 



— 1 

3- 
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Hence 



3f 
W{v° - 1) 



~2-/%Ap)^(4>) 



^WM^W+l . . (58). 



It follows from (55) that 



3z* 



1 + BIB - 2/3 



9, Q 



Wv 



2B 



Whence 



1 + /3 



+ i )WM <&*(<!>) 



1 + B (B + 2/3 



4i? 



1 + /3 



1 ) »** GW (</>) + 1 

(59). 



* * 9 9 



3 



/«; 



(.V 3 + * a ) 



2B \ ^ v 1 + /3 y 

, 1 -f if 

+ 45 (1 



PaOO<M# 



^f+f )^^)^W + 2^ 



1 . 



3a; 2 



(60). 



This is needed to express the rotation potential ^(if + z 3 ). If we add ^ to this 
we have 



\(a* + f + **) = l±^ . *±±Z l^ f ^^ 



Z 



25 " 1 

1 +5 B - 1 -f- 3/3 



45 



1 -0 



W(h>W(<I>) + ^ 



2 

1 - is 



This expression will be needed hereafter. 



. (61). 



Returning now to the formation of the expression for P/p 2 , I find 



3& 3 



1 



p* 



^~-W'-W® L 



1 + B vHB + 1 - 3/3) - (i? + 1 - £) 
2j6 > "" "T^-'yS VzW^iW 

1+ if ^(if - 1 + 3,3) ~- (5 - 1 + /3) 2 



4if 



£ 



Q,.4 



"4" 3v 



4z^ 1 + £" 

f- /3 + 1 - /3_ ' 



On considering the forms of the functions $ 2 (*>), P 2 2 '(^) ? ^ i s found that this result 
mav be written thus : 



7 9 
/i/ 



| 8 (0WW 



1 4- if 



3v«(»» - 1) ( v a - \±&) ■ 3 (1 - /3) L 2# 



1 + i/ 1 2 oi)© 3 (?) 3/3 woww , ; 
!,(,) " t_ 25 woo i_i 



Therefore, writing ^ (fj-)€ (<k) for unity, the surface density of the focaloid shell, for 
which v — v , is 



PP 



■1±B ^O*)®^) /3 WOWW , i M / w r /n 

62? p 3 (, ) + 25 i/K)" + 7?W/0«W) 



"I 



. . . -(62) 
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By means of (51), we now at once write clown the external, potential of the ellipsoid. 

It is 



I: 



i + " ^' } %M€M> + % ^WM«m + *f ' 



2/i (&, 



•IB ®., : 



. (62). 



In this expression M. denotes the mass of the ellipsoid, and the ({Ffs are merely 
coefficients determined approximately in § 10. 



In order to find the potential internally, let 



')•■■ — x- -\- y- -J- z" ; 



and, as suggested by the form of (61), let 



r. 



o 



I. + B 77 + 1 



/3 



1 ~\- B B 



1 2 /> 



-1 + 3/3 # 3 V) 
1 - /S W(x) 



WM&fit) + V 



o 



3(1 - /3f 



Then. r,j~ is a solution of Laplace's equation throughout the interior of the ellipsoid, 
and at the surface, where v = v , it is equal to x~ + !/~ + £"• 
Now consider the function 



V = — -|7rp(r- — r -) + 



"1 + 7/ ® :3 (, ) 






. 27,' © 3 f »>„) 
3/3 ©,»("„) 



fl>*(")lM/*)®*W 



®o (*o)" 



fS-^k)^ 8 w^^)^w+ S j 



• (64). 



The whole of it, excepting the term, in r : \ is a solution of Laplace's equation for 
space inside the ellipsoid, and the term in r z gives V 3 V == — 4irp> Also at the 
surface, where v = j> , this expression agrees with (63). Hence we have found the 
■potential of the ellipsoid internally. 

The potential at an internal point does not lend itself to expression in elliptic 
co-ordinates, but it may be given another form which is perhaps more convenient. 

In our present notation the well-known formula is 






1 



uC~ 



y, 






c//7 



"o 



/^-llfl) /.- d (^-.l.) /-V/V 



?)V 



■i)» 



Since $ (v) = 1, P/ (f) = ( v 



I. +/3v' 



1 



f L '(,) = („*- l>', W l (v) = v, tl 



)C 



x . -13 

integrals may be expressed in terms of the Q-functions, and we have (omitting the 
divisors (3* and E for brevity) 



VOL, CXO V JUL. -—.A. o U 



, (65), 



5 1 4. 
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In this we may substitute the expressions for ocr, y : % z 1 found above. 
It mav be worth noting that 






p7(.) t W(-) ^ *!(") ~~ v(r - i) J (- 3 - I-?)* ' 



Also P, 1 WQi'W + ^i 1 (") ©i 1 (") + #i(")&i(") + fuW€ (^) 



0. 



Tills last follows from the fact that if a, ?>, c are the axes of the ellipsoid, and 

in 



if ty denotes the function 



• n ABC 







t p7 m (proportional to our <j^ (V Q ) ), U' is a homogeneous 



function of decree — 1 in a, 6, c, and therefore 



a — ,-- 4- /> 

da 



dl 



r d x l r 



*L' . 



§ 13. Preparation for the Integration of the square of a surface harmonic 

over the 'Ellipsoid. 

If it is intended to express any function in harmonics, it is necessary to know the 
integrals over the surface of the ellipsoid of the squares of surface harmonics 
multiplied by the perpendicular on the tangent plane. 

The surface harmonic has one of the eight forms 



V = m (?) or P/ (p) ] X 






and the P-functions are expressible in terms of the P's where 



™=^(£)'>-^ 



I shall in this portion of the investigation frequently write /x = sin 0, and shall 
omit the /x or 6 or <£ in the P-, 0-, S-functions. Also I may very generally omit the 
subscript i> as elsewhere. 

If dcr denotes the element of surface of the ellipsoid, and 

1 + £\* 
1 - £/ ? 



M = Ifiviir 1 — \f[v l 



so that 1 7T M is the volume of the ellipsoid, we have, by (50) of § 11, 



pdcr 



a 






(l-/3cos2<£)* 



/1 + jS 



P 



rr 



Thei 



j,(V,r-Ar = M(l - «.f) ^=-^{^^(V;)«W+ 



where the limits of 6 are \tt to — |-7r, and of cj> are 2tt to 0« 
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It will be legitimate to develop pdar in powers of sec 2 up to any given power, 
provided (V/) 2 involves as a factor such a power of cos 2 # that the whole function to 
be integrated does not become infinite at the poles where = ± j^ir. 

I shall at present limit the developments to the square of /3. 

We know that P* is of the same form as ^*, but with the additional factor 



Suppose then that 



and let 



m 



rhen we put 



and 



y 



i\ 



cos 3 6 



n + $\\ + /3 2 H 3 = ($*)- or 



cos 3 6 

\±£ - Ti n 3 e 

l - p 



( Pf ; 



cos 2(f). 






COS 3 



/'\ • 



Now suppose that A 2 , a function independent of 0, denotes one of the four 

(CD*) 3 or (§b*) 3 cr (O) 3 or (S'j 3 
"(1 - f3 cos 2^ ' 

Then in the cases involving ^-functions and P-functions respectively, we have in 
alternative form — 



for < 



p 



\p(N?)klo- = M(l - >S)* f f jfifs- 



dud(f) . 



If it be supposed that the development in powers of sec 3 is justifiable 



7<\ 



S0 



1 + 



7 (0 + /3 3 )' 



cos 3 # 



1 



/3 + /3 3 3 /3 3 • 

cos 3 "^ a COS' 1 



n + y 8n 1 + j8~rt. 



H cos + /8 



" n o (7 - 1) 

cos 6 



-\- n, cos 



+ /3' 



•n„( 7 -i) , n„a-7) 



COS # 



COS 8 # 



J_ Uli^TL 1 ) I TT 

cos 



s0 



And F 2 has a similar form, save that y + 1 replaces y — 1, and y — | replaces 



3 



- y. 

It is clear that unless IT is divisible by cos 3 and II 3 by cos 0, [F^W and \Fod0 will 

have infinite elements at the poles, and the development is not legitimate. 

cos 5 6 d i+s 
Since P* = - — - -y- (/x 2 — - 1)*, it follows that the power of cos # by which P 5 is 

^j . % \ CL Lb 

divisible increases as s increases. 

3 u 2 



f>i.<; 
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Let us consider the case of s = 2. 
Then 

ii + /3n, + /3% = [P* - ^ P - p qi v* + WJ 

= (P 3 ) 2 - 2y8 (g PP 2 + 9iP 3 P 4 ) 

+ & [2 ?fl P*P« + (%P) 2 + (ftP*) 8 + 2 M4 PP* 

(or the same with accented qs for the other case). 

From this it is clear that IT is divisible by cos 4 and II x by cos 2 (9, and the method 
of development is legitimate when s = 2, but it is not so when s = and s = 1. 

The investigation then separates into the general case, and the cases s = 0, s = 



§ 14. Integration in the general case. 



We have 



and 

(if 



r = p __ ^ _ 2 p - ^ _ ^ f 2 P . + * + ^ , ( p - i + ^ + 4 p + < , 

= (Pf - 2/3 ( ? , .... , PP - * + <j ! + 2 PP + 2 ) + 2£%_ 4 FF-* + ? i + ( PP' M ) 
+ £ 2 [ (g. _ 2 P' " 2 ) 2 + (?, + 2 F + 2 ) 2 + 2?, _ o ? , + , F - 3 F ' 2 ; 



Also 



.. / cos 3 



■F 






-0 



sin 3 



has the same form with accented c/h, so that it will be 



merely necessary to accent the r/s to obtain the second case. 
We have then 



n = (F) 3 , f 



2 (g,, _ 2 FF - * + y, , , FP ! 2 ) , 



n, = 2 (? , _ 4 FF - * + g , + 4 FF + *) + ( ? . _ 2 F - 2 ) 2 + (?, + 2 F + 2 ) 2 

+ 2q. - 2 q. + 2 P' 

Then since cos 6d6 = dp, 



s - 2 T)s •!- 2 



I 



in 



FdO 



«' — }. 'IT 



\vfdii + j3 f + ' ^- 1 ^-^dp - 2$ [ + \q,^ i PT'-- 2 + 7, , ., FP* * -) dfi 



-1 



1 — ylF 



p ] (7 - 1) (P) s , , m F (I - 7) (E)= , 



9 oof + 1 (7 - 1)(?.,- 3 P S P S "- + ?. ! + 3 EP' H 2 ) , 

2 H-i" " " i" : F 



^t 



+ i 



+ 2 [2?^ PI 



> - 



. __1 



+ 4 



+ 2q, ., 4 FP" + 4 + (r/. _ 3 F - 2 ) 2 + (?. , 2 P + 2 ) 2 + 2r /s _ ,qr. + 2 F -^ 3 F ! 3 ] ^ 

. (6G). 

F o d0 has the same form, but with accented q'&, and with y + 1 replacing 
y — 1 , and y — -| replacing | — y. 
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5 I 7 



It is now necessary to evaluate the several definite integrals involved in this 



expression. 

It is well known that 



(P/(F))^ 



2 {_+_*} 
2i+ \ ' T— ,s! 



It is easy to see that it is possible to express P**-* in the form 

P; + ** - AP, + BP/„. 2 + OP;L , + ..., 

where A, B, C . . . do not involve /x. 

The value of A may be found by considering only the highest power of /x on 
side of the identity. 



Now 



1 



"3 * + 2/t 



(i -^fc+A/rf v + i, + 2 * 



2* •? J 



,rf/V 



(F- I)'. 



(-)*•■ 



- & 



2i! 



2» . i ! i - § - 2/: ! 



i /^ "i • • • > 



and 



2i\ 



V* = ( — ^ - uj 4- 



2 % .%>■% — s\ 



n\ 



I herefore 



xx — ■ ( ~~~ ) T" 



I — Si 

7 . ^"ZT^T 



Then, since the integral of the product of two P's of different orders vanishes, we 
have 



+ W +3 *^ - (-)* . 



% — - s 



o __ o __ 91- ! 



'V) 2 ^=(-)*57 



9, 



% 4- s 



» t 



+ 1 * — s — 2& ! 



P*P* + 2& 



We will next consider r.-— — ~ei/x, where & is not zero. 

J I — fl' 

The differential equation gives 



d 

dfji 

a 

dfx 



(1 - f) 



d¥* + 2/; " 



^ 



+ i(i + 1)F 



•!- 2/,; 



(s 4- 2A) 3 p .,. 2/ , 
I — /x* 



dP* 



—^)^ +*(*+l)I 



3* 






I: 



/^ 



0. 







Multiply the first of these by P s and the second by P" ] 2/ and subtract, a 
have 



we 



fc (> -f k) P'P« + 2A 



F 



+ 2X- 






(1 - /**) 



tZl 3 *' 



Ptf 



f ¥ . 



(l - ^) 



).s' + 2k" 



dfM 







18 
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Mi 



Therefore* 



'+ 1 T>s]M + 2k 

4k (s + k) | ,—;., dfi = (1 - /^) 



J _ 1 1 



Again since by (11) 



r 



\M + 2k 



it follows that 



I 



> + 2& 



ffyt 



0. 



Pi 






Vc 



between limits i 1, 



A p, -,- 2/, -,- 2 + B p s+ 2/, + p,4-2/ ( :.-2 ? 



>+1 P*P* + 2/t ' 



d/x = 0, unless & = or 



It remains to find the integrals of --- --— , ■r—-^~irz 9 and ----- --- 

° J. — fi j (1 •— //*) : * (1 — jj?y 

We have Been in (11) (transformed to accord with our present definition of P*) 

that 



P* 

1 — /X" 



XXClivv 






4-9 (s -f 1) 



D* + 2 






I 

9 



*-i('i -f 1) 






1 



amjifl'iWi 3 



P* 4- ^'J-LSli 8 r j p? 

^ is (8 - 1) 



V\2 



(P') 



l 



/*" 



p Cv/X 



1 



4-9 (s -f 1) J 



pp + 3 4t + 1 



^(t 4- 1) 



«**i 



+ i 



(P) 3 C?/x 









i 



4?fe + 1) 






^ 



i (t + 1) 

O -1 






{i, s} {i, s — 1 } f 

\ / 

(p) 2 



P*P* "" 2 tZjU- 



1 



p 



V ^/^ 



+ 



{i4R* - 1} rvp*-* 

4s (s — 1.) J I — /a :: 



a 



f LAr , 



■ p^p* + 2 

(1 ~Z ^iji r/ ^ 



1 



(p + *)* 



4<? (.9 + 1 ) J 1 — fM 



- da + 4 



*(* + 1) , -1 fpp + 2 , 

*~ — .1 Jl ™ /T ' 



[i, *j- J*, S 



H % s l - $*■ 



11 fP" 2 P 



4s ($ — 1) 



O ( f fA ,- 



M 



The first of these involves integrals already determined : on introducing them on 



the right and reducing we find the result to be ~ . -~' 

° 8 'I — S I 



The first and last terms of the second integral vanish, and the integral is clearly 






"i{i ~f 1) 



S' 



1 



+ 1 



1 i + 6' ! 
s % — s I 

The second and third terms of the third integral vanish, and the whole is clearlv 

1 i + s -f 2 ! 



r*r 



4s (s + !)(« + 2) *_ 8 -2r 



* I owe this method of finding these last two integrals to Mr. Hobson 
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Collecting; results we have 



+ 1 

1 

+ 1 



(P) 2 d/* 



■i + * • 



% + 1' i 






7' -f- eS ! 



2* + 1 >■ 



- ] 1 — /i- 



0. 



,M-1 (J*);! 



- I 1 - jLt 
+ 1 (p,)3 



- ~ rf/X = 



1 i + s\ 



8 % 



S : 



f 



1 (1 - fj) 

1 1 p*p* + 2 



2\2 Ct/X -g 



(1 - /* 3 ) 2 

+ 1 ^>^i> 8 ""* 2 



d/x 



"i(i -j- 1) 



+ 1 



2\i ^ 



Then by means of (66) and (67) 



2k I 



1 i + a ! 



si 



i ~f $ -f 



4s (s +l)(s+2)' i - 5 

1 i + sj 

i(s Z 1)^3*2) * i ZZ7! 



2_1 

9 ! 



/' 



7 



• • ( 6 ?) 



/<>/# 



2^ -f 1 i — s! w 7 .s / — *! 



+ 2 /3 &„ 2o . 






2i -f 1 



/ -r 6* 



7, 



9 ! 



S 



~r (/* + 2 • <> • 



+ j8*(y- 1). L . ^4| + ^(f -y)| 



s i 



21 -(- I 
»_(* + 1) 



i -f- a i 






+ 1 



1 i_4- * 

-s " i — 8 



. --Mi 

"*" 2i + 1 



% - 4 



% \ s 



i — s ! 



4 ! i + s ! , \2,i + s 



2! 



£ 



6' 



4! 



Therefore 



Fd0 



+ £. 



+ fe + 2) 3 T 



i + s +• 2 ! 



s 



7m + </--2^ 



s + 2! 



+ 2 



I 



S 



i + 8 



2* + 1 ' i 



S 



*"l / • 






,: + ifff7-T) + «*.-(*-•»• 



_2 

"9 ! 



5 



1) 



Wq 



S - 4 



(i + «) (■£ + s — 1.) (i + s — 2) (* + s 

+ 2/8y + i (i - s) (i — s - I ) (i - s 
, /3 3 (q. - 2)- 

i~ i; 



3) 
2) (i 



6' 



3) 



(i-8+lj(i - S + 2)(* + s)(i + s - 1) 



+ P*(q t + i )*{i + s + 1)(» + « + 2)(i - *) (*' - s - 1) 






I (' + » ! 



1 r 



I - .s ! 



(1 + /8) (y - I ) + (3 (| - y) . i p| -^ + 1.1 }. . (68). 
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,1 



Also iF^dd has a similar form with accented qs y y + 1 for y — I, and y — \ for 
When we substitute for y its value I — cos 2</> ? and write as bo fore 



v 



■/ (V + 1 ) 



•S" 



the last tern:.! in FAQ becomes 



1 i -f H ! 
+ /3 . . ;. ~j {i/8(S + 1) — eos2</>| 1 



1 

9 



S{t - i); 



, (68), 



Also the last term in FAQ becomes 



2 l 



+ £ . : . v- - ^ {2[1 + \p(t + y)] - cos-2</> I + J-/8(S + 3)]}. . (68). 



* 'I 



8 ; 



But it will appear later that we only need the parts of these terms which involve 
cos 2(/> developed as far as the first power of /3 ; hence in both cases we may write the 
latter term inside { } simply as — cos 2c/>» 

Our general formulae for the q coefficients apply for all values of s down to s = 3, 
inclusive, although the result for «v = 3 needs proper interpretation. Hence the 
present result applies down to $ = 3, inclusive. 

I have just re-defined S, and I remind the reader that 



r 






Then if in (68) we introduce for the </s their values, we ihnl that the coefficient of 
the term in /3 is 



'^-z 



( ,:+7)"(i"+i-l) + tf* + * (*-*)(* 



tS 



ft* i A«rt * / a 



The coefficient of the term in /3 2 is 



64 



( i— $ + 1) ( i— s + 2) (i — s + 3) (i — a ~h i) . (I — s) (Z — ^ — 1) (/' — .s — 2) (i — -s — 3) 

"($ - l)(s - 2) ~*~ "^ ~ " "~~ — <v^H^ 

(i + *)(*_ + $ — 1) (j ~" $ + 1) ('/ - a -f 2) + s + 1) (i + s + 2) (i - 6-) (I - s - 1 ) 

(i — -s -f !)('& — g + 2)(i — -s)(i - * — 1)| 



If this be reduced by a process similar to that employed in § 10, we find 
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F,<W 



„ , t , | I *■"*"■"■ i^fji^i ~~" I ) ~p 32/5"' oJ^" ~™*" 02* ~y~ Q ~~" tS w ( 2/ *~p Z^ *""■ 1 j ~j~ JI j 



1 



— j— jo 



1 i + s I c , 



s % 



- {£0(2 + 1) - cos 2<j)} . (69). 



We know that P y is derivable from ||J* by multiplication by 1/C/, and we have 
found in (33), § 9, 



1 



;7i = 1 -f j3(t + 1) + ^[3S 2 + 8S + 5 + 5 2 S- - T 



(CA) 



Hence multiplying (69) by 777-- we have 



F o d0 



2 i ~f §'! 

2?: +1 ?; -7i 



(Ctf 



(1 + £0(2 + 3) + a ] #[- 2 3 + 26S + 42 



+ ^(82* -. 2S + 1) - 2T]} + £• 7- ;.-_--; {£0(2 + 1) 



1 i + s ! 

e9 1, — S ! 



cos2</>} 



. (70). 



I have also obtained this result by direct development. It may be thought 
surprising that the last term is now the same in both formulae, notwithstanding the 
difference in the earlier stages, but if the reader will go through the analysis he will 
see how this has been brought about. The formulae (69) and (70) also hold true when 
s = 3 (as I have verified), notwithstanding the fact that P 8 is not to be derived from 
P 3 by a factor. 



The next step is the integration with respect to c/>. 
We have 



< 






^ s 4> + ft [p» -a j °°* (* - 2) </> + p., h A ^ (s + 2) $ 



cos 

si 



+ P 2 l;p^^ T (s - 4)r/> + ;>,, 4 |^ S (5 + 4)r/>; 



Therefore 



1 m* 



- 4" ± I cos 2*^ + /3 [(jp„_ » + p, 1 2) cos 2<j& ± 2>*~2 cos 2 (s — !)(/> 

± P, + 2 cos 2 (.? + 1 )$ } + /3 2 [|(p^ 2 ) 2 + £(#h-2) 8 + (p,-4 + p, 11 
+ p^2i>, + 2 )cos4(/> ±p,- a ^ + 2 cos2.90 ± (p.^4 + |(p^o) 2 )cos2(.s s 
± ( T>* 1 4 + M P, + 2~) 2 ) cos 2 (* + 2) <£] . 

Also i '-—zr — —-■ have the same forms with accented p\s. 
.1 — p cos 2cj) L 



"— ■* *j / (7.) 



Accordingly, with unaccented f/s, we have to multiply this expression by 

VOL, CXOVII.— A, 



3 X 
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(1 — fi cos 2<f>)~* 9 and with accented p's we multiply by (1 — /3cos 2c/>)% and we shall 
then have the functions denoted above by K% 

The function K 2 has to be multiplied by a function of the form A -f B /3 cos 2(/>, 
and integrated from <jf> = 2tt to 0. It follows that the only terms in K 2 which will 
not vanish are those independent of </> and those in cos 2</> ; moreover, the latter terms 
are only required as far as the first power of /3. 



N 



ow 



(1 - £cos2<£)-* = 1 + \fi cos 2(j> + fV^(l + cos4<£) 5 
( 1 — £ cos 20)* =1-^8 cos H — ^^ (1 + cos 4c/)) . 



Then omitting terms which will vanish on integration 

(1 - ,8 cos 2^j 4 ~" 

(C/ or S/) 2 



1 — /3 cos 2<f>f 



\ {i + i8--[(p.. 2 ) j + ( Ps+2 ) 2 + j P ,_ 2 + i/>, +H + tv; 

+ fi( P* - 2 + .P, + 2 + i) cos 2<£ , 
+ i8(pV~2 + pVi -2 — i) cos 2^. 



1 
i () 



However, the latter formula is not needed except for verification, because it will be 

1 
derivable from the former by multiplication by 7 ^-v, . 

Now if we substitute for the p's their values as given in (27), § 8, we find 



(«!V_or &£)* 

(1 -£ COS 20)* 



M 1 + A /3 2 [S 2 + 4S + + s 2 (S 2 - 2S + I)] } + \l3(t + I ) cos 2c/) . 



And multiplying by j^~^ or 1 



jift z (% + 3), or developing directly 



(C^orjS/) 2 
(1 -/3cos"20)* 



M 1 + ^yS 2 [2* - 6 + ^(S 2 - 2S + 1)] } + }fi(t + 1) cos 20 . 



These represent the K 2 of our integrals. 



Then 



P[\% 



. [ c- 



8\2 






d 



(J 



M 



(1 - 0f f |>j v^^'-^r; fW# , 

v ; J J J n. — /3cos2rf>v ' 



2wM (1 - /3> ! i + * ! , ., 



2/ 4 1 



7, 



1/3(2 - 1) 



+ ^?)8 3 [3S 2 -~ 62 + 6 - s 2 (2 2 + 






1) + 2T]} 



X {1 + A/3~[S 3 + 42 + 6 + * 2 (2 2 - 22 + 1)]} 



+ i 7rM (l- i 8)' 



^.l.*^Li; (s+1) „^.i..*+l| (s + 1) 



,s s 7 -— S 



.s -i — .<? ; 
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2ttM i+j ! , 
■li + 1 i - s'l \ 



■S 1 



2ttM % + a i 



2i -f- 1 i - .s ! s 



+ iViS 2 [22- - 2 + G - (r(2% - 1) + Tj ) 

■ i/SV. + f L/3 2 [22 2 + 32 - 6- (22 - 1 ) + T] } 

.... (71). 



The following results may be obtained either by direct development, or by 

multiplication by either or both the factors -r-rz; and —-— -. . The former converts 
1 J {(iff (D/) 

!P into P, the latter C or j$ into C or S. 



c 



A3 



P[^i\$\)d<r 



27rM i + ^ ! 



2i + 1 i — 5 ! c * 2H " 



{1 



+ -± € p*[2t* + 2 — 6— * 2 (22 — 1) + T 






/ 



p ( p/ 



°0; 



da 






i + pes + 2) 



+ 1T#[HS + 10 + 5- (22- - 22 + 1) - T 



f 



3' ? 



r / 



plP 



■< 






* * « 



(1+1/3(2+2) 

+ ^r/8- '[9S + 4 + s-(22 2 



25^ + 1) "*"*■ TJ j- 



y 



(7 I ). 



§ 15. Integration in the ease of s = 2. 



Although the development in powers of see 2 is still legitimate in this ease, yet 
the formula? found in the last section, fail because T contains s 2 — 4 in the denomi- 
nator. Moreover since :J$ 2 is not convertible into P 2 by a factor each case must be 
considered separately. 

We now have ^_ 4 =0, ^-.,^=0, and therefore from (68) 



F x dd = 



2 i + 2 ! 



2-t + 1 * — 2 ! 



i 1 1 + 2i8 



2o 



_(* + l)(i + 2) 



+ /S 3 



+ <lS - 2)(» - 3) 



(So> 



2 2o (* - 2) (* - 3) (» - -1) (,■ - 5) + ( . ~~ . ( . + ( . + 2) 



+ W (* + 3) (* + -1) (i - 2) (i - 3) + 2 Ml , ^§| T f 

V mm* • 



/'W0 is equal to the same with accented </\s, and the last term equal to 



i + 2 ! 
i.^g_...__rj 1 2 + |-/3(2 + 9) — cos 2gf> 



O X. ^u 



«J «w "X 
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We now have 



</o = 



-J-{i, 1 1 {i, 2}, <7 4 = ^4, <y G = ysVg? associated with a cosine function. 







<?'i 



x ] 2 5 </ fi = 5 } 2? associated with a sine function. 



It is well to note that these values are given by the general formula, because this 
consideration shows that much of the previous reductions is still applicable. 



Effecting the reductions I find 



F,i 



9. 



i -f 2 1 



(i - kP (t - i) + ( hF l.¥^ - ¥* + 20 



2i +1 i — 2 ! 

+ 2 (2i + 1) (35 - 1 )]} + h^lL'll i^ (* + X ) ~ cos 2( ^i 

2<f+ T i~ 



i + 2 



-■ {1 - |/3 (V. - 1) + 2^^(19^ - 1302 + 80)} 

?' -f- 2 ! 
+ -J/3 y— ; ^i(i/3 (55 + 1) — cos 2<j>). 



V 



This integral will be associated with C/ and C,- 2 , aud in the present case 
= 3 I [I- -f- 1)« 
In the same way 



,F o d0 = 



9 -/ 4- 9 ! 

2i + li-2! li ^ 2/5 



;S + 8 « (2*' + l)j + fep* l^i-t : + ¥S + 9 
6(2 i + 1)(S + 5)]} + i/3 ^4r i 2 + M- + 9 ) - <*>s 2 <^< 



9 / _L 9 f 



*•*' ' *| J, ' -■■'■■** ^ „ 



5 -; { 1 + P (S + 3) + - 3 1 - s /3° (25S 2 + 1862 + 368)} 



z' -4- 2 ' 
4- -Lfi ----- < f 



$p(% — o) — cos 20 j. 



This wi.ll be associated with <S 7 2 and Sf» 



These factors were 



Now turning to the cosine and sine functions, we .find that they must be treated 
apart, but the integral involving Of may be derived from that in C 2 by the factor 

2 (cos) ; and similarly Sf from cS 2 by the factor —(sin) 

evaluated in (36), § 9. 

We now have p g ^_ [t = 0, p\_ Li = ; also for the sine function ]),__., = p, } = 0. 
Then 



J 

n 



mr 



2' + |- cos 4<£ -j- )8 [(2jp + jj|.) cos 2f/> -f- ju, t co» 6<^> 

+ F[(PoY + UP-tY + (P« + i'uPi.) cos k/> + p^coa 4</> 

+ (?>« + K/'i)') cos ^ 



(Sf)" 



2 



2 



cos i<£ + /3 [ p^ cos 2(f> — p 4 , cos 6<£ 



+ ^Li(^) 3 + lh cos 4(/> - (^, + 4(p + ) 3 ) cos 8<£J< 
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Then as far as material 



(®/) 2 
(1 -/3 cosily 

(1. — /3cos2</>) ! 



Now 



(1 



(«/)* 



- /3 cos 2<£) i 

(I — /3 cob 2^y 



iO- + F\'{/hY + 2U) 3 + Po + -^ + -A-]} 

+ fi(p h + 2p (J + I) cos 2(/> , 
Hi + P Kpd* + £/>.,. + A] } + /3(p 4 + i) cos 2</> . 



ri — 1 s -1 9 '. 



24 p 5 Oj and 



■Hi + ^^'(l^' - 8S + 22)} + i-p(5$ + J)cos2</>, 
H' 1 + \hP*@' 2 - 81 + 18)} — &8(2 - 3) cos 2c/.. 



We now multiply these by i^eW and i^(i0 respectively, and the last terms dis- 

appear as before. I remark that the disappearance of the terms which do not involve 
the factor l/(2i + 1) affords an excellent test of the correctness of the laborious 
reductions throughout all this part of the work. 
Then we have 



p('PfW) 2 '^a* ; 



9, 



IT. 



i + 2l 



2i+l i—2\ 



2ttM 1 + 21 



2i+l t~~2! ( 



j(l -£)*[] ~|.)8(S-1) + - 2 - W 3 (19S 2 - 130^ + 80); 

X [1+^^(19^—85+22); 



(72). 



If we multiply this by 



D 



~(cos) 



or 1 — y£-/3~(55 + 7), we obtain the result when 



G'i replaces (&? ; the only change is that the last term inside { ] now becomes 
Again 



+ 1 '(' 



x|;").+- < r4-/8-(S--8i4-l«X 



_ 'I -j" X 'I — Jj I 



(72). 



If we multiply this by _'"'., (sin) or 1 + o/S 3 (S + 5), we obtain the result 

when St replaces jg/ 3 I the only change is that the last term inside -j j- now becomes 
+ Y l^ (295- + i 065 + 136). 

This terminates the integrals, which can be completely determined by this method 
of developing in powers of sec'- 0. 
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§ 16. Portion of the Integration in the case of s 



The preceding method may be used for finding the four integrals 



and 



H<cK< p 'M p ' V ; %; 



sin 3 # 



D™ 



w. 



a 



There will then remain four integrals of the type p (S/)~ ( P*)' 2 d<r to evaluate. 

The first pair of our integrals are dearly to be treated by putting II = 0, 



q s _, = q m] = 0, q s „, v = g_ 3 = 0, and then determining ¥ x d6. The condition for the 
second pair only differs in the accentuation of the qs which vanish, and in the use 
of \\d0. 



The vanishing of II makes 



F-, ^IL cos 6+fi 1 ( n ^-.i } +n, cos 0) 

1/1 \ (308 ■* / 

F^/Slii cos #+/3 3 f 11 - 1 .^,^ +1L cos 6^ . 



COS 



In the first of these 



and in. the second the form is the same with accented (/s. 



I )l P> 
Also since It — - a 

1 — yL6 " 



fii,rf0 

0, we nave —■—==(). 

J cos 6 



H 



ence 



F 1 ^=-2/^ 3 |FF»^+2/3M[( 9D Pa«+Mr i3 )=- , (l M )->//*, 



2 Vy -j~ 1 i 

27+1 ^1~! 



;2/8j 3 ((' — :i.)(* — 2)-+- 20% ( 4 - 1 ) (i — 2 ) (i - 8 ) (i - 4 ) 



+iS a (g 3 ) 3 (*+2)(i+«)(*-rj(*'-a): 



aire 



1 



i\ 2 dO is the same with accented q } s. 



lb is only necessary to pursue the cases p (j|/)fF,<:/cr m.k.1 p (C; r )' 2 Fo^cr ? since tire 

other pair of integrals may be determined by means of multiplication by the 
appropriate factors, determined in § 9, 
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It 
/ 



Now for p.; 1 $ ? -\ associated with F 



1.3 



2s=i 3 &[1-t&W+1)], ft=rk; 



and for P/C- 1 , associated with F 2 , 



f/3 / =i(r[ l +i L ^( i + 1 )l ?r/=T-o-8- 



m 



Therefore 



F,eZ0: 



F„cZ0 



_2_ -i + lj 

2?;+ 1 iT-^T"! 



2 i-fl! 



2?;+i *-i i 



W-l)(t-2) + iki8 B ri(»-l)(i-2)(t 
+^(i+2)(*+3)(t-l)(i-2)-(»- 

^(*-l)(»~2)+ T b/8 s [4(t , -l)(*-2)(t- 
+ -f(i+2)(*+3)(i— l)(t-2) + 3(i- 



■3)(*-4) 
l)(»'-2)t'(» + l)]} 



3 



)(t-4) 

l)(z-2)i(t + l)]}. 



In the present case we cannot use 2 as an. abridgement, since it is infinite ; I 
therefore now write 



Effecting; the reductions we have 



fa< u =x 



2 i + 11 



2i + l i — 1 ! 



:{P(y— 2(2»+l)+4) + T bi8 a [-«V 4 +Vy + 30-1fi(2t+l) 



h 



F ^=a.- +1 t 



2 t^M/8(.y-2(2i+l) + 4) 



+xh/8W+ L ! s /+ 174-16 (2»+ !)(,/+ 5); 



The former of these is associated with <£, the lattei* with (£. 



In the cosine and sine functions we have 



(« 



o 



|- cos 20 + /8 (Ps cos 2 — p : . cos 40) 

+ yS"[i(Ps) 3 + .p s cos40 



and 



(P-O + I (Ps) 2 ) cos HI 



(C; 1 ) 3 = i + I cos 20 + y8 ( p 3 cos 20 + Ps cos 40) 

+ i8 a [i(p 3 ) s + R- cos 40 + (p- + £(p 8 ) s )cos6£ 

As far as material, we then have 



(a/)' 2 



(1. - /3cos20) i 



= *{1 



i{l _x/3 + /3*[( Pr5 r + i P;i 4-VV]} 

+ (~ \ ■ + i0 + « cos 20, 
- 1/8 + /3-[(p ?) )- + iPs + fVIKl ~ (1 ~ ¥pa ~ i^) cos 20}; 
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(©«')* 



(1. - p cos 2£) 4 



|{1 + i/3 + /3~[(.p ;! )" + hlh + rd} + (4 + ££ + #p 3 ) cos 2f. 







{1 + i/3 + /3~[(p 3 ) 3 + £#, + A]} {1 + (1 + 2/3p s + i/3)cos2<£} 



<j 



Now both for sines and cosines, to the order necessary for our present purpose, 
— ~ is {h 2}. Therefore, introducing / for i(i + l), 



(&; 1 ) 3 



(1-/3 cos 20)* 
(1.-/3 cos 20)* 



1 {1 ~i^+^k^(.f- 12/+ 68) J- { 1. -[1 +A0/--&8] cos 20} 



1{ 1+1/3+ s!jf)8 2 (/ s - I2/+G8)} { I. +| [1 -J / Q/+i / Q] cos 20 



r* • ("1). 



? j 



Observe that|Fjd!# and F.,d# do not involve cos 20, and are of the first order in /3, 
II once, as far as material for \he premnl portion of the work, 



JS2 1 \ -7 



(fe 1 ) 



s 9AU 9 V 1 - 4P/: 



( 1.-/3 cos 2<f>) 



Also, to the first order, from (37), 



~ (sin.) 



tew *' 



j TjT /3 « 



l^iZ — li'-i J-.im 

(I -£ cos 20)* 2V A T"4P; 



IX- 1 



(cos) 



J—' I 



(73). 



Therefore as far as necessary 



.(--$* -.l/i _l_ kg) iOH 

(1-/3 cos 20)* 2 V ■ " ^ 4P/ ( I . - /3 cos 20)^ 



Hence 



i- ( I " ;{,/3). 



?> (*;T CGV) 2 



(W= I?* ; + M (I -/8) 1 (1 H/3) {^ (/-2 (2/+ l)+4) 



"7 + I i 



+ x a s/3 2 [ — <i/ ? + W+ 3 ° — 1 6 (2* + 1 ) ] 






(74). 



For S/ we have only to replace the factor I — -|-j8 by 1 -++/3, and find 



2>(S ! ) ) 2 [i/) 2 ~(I.V) 2 ](icr 



^HtT^O'+^-tA^O" 3 






32; — 84)} 



• e 



(74 
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Again, omitting intermediate steps, 



pwyUpw-Iw V^h-^ e J 



p(on 



d 



(T 



2ttM i+ll 
2i+l i-ll 



_ 7rM i±lJ [ l /8+ i /83(2i+7)] 



(P.l)2_(p 



\ 



V 



\ O' 



[« - sin 3 0Y 
cos 3 



rJ 



a 



^ mmj+i)+fh0 ^ +ie , +im 



TT 



M-. 



I -f* 1 : 



•~[|/3+i^(2y+i)] . . . . 



■ • . (74). 



» » » 



(74). 



§ 17. Portion of the Inter/ration in the ease 0/5=0. 

We are to find \p\ *\L [(P/) 3 — (P«) 3 J ^ leaving two integrals of the type 

lp(C 2 ) 3 (P«)^cr to be determined subsequently* 

It is only necessary to consider <&, since the other case is determinable from it by 

1 

multiplication to "pprri, as found in (88) of § 9. 

Following the procedure of the case where 5=1, we have 



Y^W 



2fa PP-^ + 2/P [? 4 PP fc +i(ft) 8 (P 3 ) 8 ]^ 



2 



2i+l 



\2fai(i-l) + 2Pqti(i-l) (t-2)(*-3) +)8«fe) 2 (t+l)(* + 2)i(t-l)3. 



licit ©J.I1UO c^n — ^j fjf| ( — l"2 8> 

Now • (C) 3 = 1 + 20p 2 cos 2^+0* [(2^+1.(^)2) cos 4</,+* (p 2 ) 3 ], 



anci as tar as material 



(1— /3 cos 2$)* 



= i+^ 2 [i^+i(p«) a +A]+i8a+2^) ^ 2 ^ 

= U+Ai8 2 -2 --4; + 6)}{l-^(y-l) cos 2c/>}., 



(75) ; 



since 

VOL. 



jP2=— ^(^+1-) 



Cv'TT.—i. 



.1 ./ 
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At present we only require this to the first power of /3, and since Fj<i0 does not 

contain cos 2<£ ? the expression (75) as far as at present needed is simply unity. 
Again, by (38) of § 9, .:-. 

1 l+i/3 3 0~3), 



(I);) 3 



therefore by multiplication 

(CO* 



(1-/3 cos 2cf>y 



= U+AW-6)} {1-P0-I)cos2^} . 



• (76). 



This is also unity to the order at present needed. 



J! L tJX.Lv.' O 



V 



[{o^iW-iPiYWr-- 



2~-M(l-i8)*Ui8[i+l-(2.-+i)] 

+^)8 2 [5/+14;-+12-4(2i+l)0'+3)]} 



^rM^ + i^O-l)} 



( 77 )< 



§ 18. Preparation for the Integrations when s= I and 0. 
We have now to evaluate the three integrals 



Ju 



M 



p; \ 3 /1+/3 



2> 



COS0/ \l — /3 



sin 3 #) c?<x, 



N = fp(CiP ? ;) 2 ^ 



/ 



and from these to determine three others when S, C replace <%, <£♦ 
We have 



ji9(i< 



, jjf(l- fl)*(ffi ] fcos^+ ggg' 



rfftty [ (1 -£-jj sin s 0)* J 1 (1 -/3 cos 20) 






• (78) : 



pl£S- B m*6)d<r r j/f/l ovxA+ySN*,, 1-/3. aflV 



> < 



C0S 2^ + ^|^) 



(1— yS cos 2<£) 



.* 



It is the second factor which alone involves ^>, and as I shall now first integrate 
with respect to <j>, the first factor may be dropped for the moment, and the second 
factor multiplied by the squares of the cosine or sine functions. Since the integration 
is from <f> = 2tt to 0, those terms which vanish on integration may be dropped. 
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For brevity write 



70=^-0'*- 4?"+ 6), 
/ 1 = <rb0' 3 -12/+68). 



Then we have seen in (73) and (75) that 



(%ij 



(1-/3 cos 20)* 

(gP, 1 ) 2 _ 
"(1-/3 cos 20)* 

(®s)» 



(1— /3 cos 20) 



4 



1(1-^+^) ]l-[l+lyS(y~4)]cos 20}, 
= l(l+|^+^./ 1 ){14-[l- 8 ^(y-4)]cos20f, 
= (l+j8Vo) {1-1/8 (/-I) cos 20]}. 



lhereiore 



-\ 



>2tt 



rt _ a A , /5(1 — C0g2$) 

^i±=BLl . (W*=*(l-i0+i8Vi) {cos^+| y 8+ 1 i 6 ^(i+20)} 

I (l—p cos 2$)* 



•2tt 



• • 



. (C i 1 )^0=7r(l+i^+/8Vi) {cos^+P+i 1 6 ~y8 3 0'+4)}, 



P W 3 ^=2ir(l+/8 3 y ){cos s tf+i8+ii8 2 (y+3)} 



. (78). 



Now pick out the parts of pdcr and of these integrals (78) which are independent 
of 9, and write 

F^M (1 -/3f (i=|)* (1 -i/S+TO 



A.lso write 



G=^M (1 -J8)*(£|)* (1 +ii8+/8»y 1 ) 

= ^{1+1/8+^- /8»0'-12; + 68)}, 
H=2^M (1 -iS)* (f=|)*(l +)8 B yo) 

= 2ttM {i-|j8+AW-4;'+84)} 



.9=1(1+^"+^), 



/>=l+i#+f/3, 



fv JL f\, ' 



1-/3 



l-f/3 



3 Y 2 



/ 



• * * • 



• * 



;, so that /<*= 2/3+ 2/3^-4- , . . 



(79). 



(80). 
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Lastly, in accordance with the usual notation for elliptic integrals, write 



1 — B 

~= I ---——• sm- c/= L — k - sin- (J 



#**?»»* 



(80). 



m 



Then we have 



L=F -^^-(P/ 1 )^^ 



;-7T 



M=G (cos 3 0+/fy).A.(P/yW 



* * * 



• (81). 



— trr 



N=Hi -i ' cos ^^(p ; yw 



$7T 



The next step Is to express the squares of the P's In a series of powers of cos 2 6. 



When P/(^) = ( — ] (/x s — ■ 1 )% it is known that 



P,-(/*/*'+(l-/t B Wl -/*'*)* cos <£) = P,-(/*)P,-(/i') +2*i , T-^Pf , (/*) P/V)cos a<£. 

5=1^ + ^1 



?,~~s! 



By putting /x^/x' we see that 2— — ; (P/(/x))* 2 is the coefficient of cos s<j> in the 
expansion of P,-(l — (1— ft 2 )2sin 2 ^). By Taylor's theorem this last is equal to 



T- i 



2 (~) r(1 -:r( 2 ^ 8 W("£ W* /^ l ) • 



r=0 



d/A* 



Now 



d V 



) P. 



1 /rfN^ 



/^/ x * 2*i! Vfyj ^~ 1 ^~2 r .Tl 



1 i + r ! 



% i — r 



1 + terms involving powers of /x 3 — 1 



Also 



9> 



1 i 4- r ! , 

, r -~;, , when /x=.L 



Z' . r i i-—r\ 



m-.Srll- 



snr r |^ 



/^V-l — (j-^V-lVSr 



2^/ — 1 



< = r 






-MflJ / « 



£ = 



2 2 '.2r-*U; 






On putting r — t = s 9 we see that the coefficient of cos s<f> in sin 2 '' •§</> is 



(-) a 2r ! _ 

2 2r ~ 1 * r-~s\ r-\-§ 
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Hence we have # 



i- s I 



2-~;(P/(fx))^2S (-) 



r+s 



2r! 



Ur! 



i, + s\ 



T~S 



2 2r r— s\ (V !) 2 r + s\ i — r\ 



ccm r . . . (82). 



Now suppose 



Then clearly 



Therefore 



(P i 1 ) 2 = S y 2 r~2 cos 2r ft 



\-L j} 



S «2 r COS 2r ft 








ya- 



, , r+J _ w. , i + r ! i + 1 ! 

V ~~ / 2~' r- 1 ! (r !) 3 r + 1 ! ' T- r ! ' i—il ' 



( 



a Zr — ( ) ( 



2r I i -f r ! 



2 3r (r!) 4 i-r 



L r ^' f^ cos 3 '' 4 - 3 + /3/cos 3/ ' ia 
-jz—% y 2r ^ 2 I — — : dV. 



■It 



A 



- — 2y 2r „2 



rJT 



(cos 2r e+/3gcos 2r ~ 2 6)Ad0, j 



N" 



r=l 



r=i 



— $TT 



— 2 a 2r i 



* cos° r +2 + @h co$? r 



d9 



t t 



(82), 



(83), 



The evaluation of these integrals depends on two integrals only, namely, 



'cos 2 * 



A 



dd 



and 



cos 2 *#. Add, and these will be considered in the next section. 



— r— d0 and cos 2n 6 Add. 

I will denote these integrals D and E respectively, and I propose to find their 
values in series proceeding by powers of k*. 

The usual notation is adopted where U(x) is such a function that it is equal to 
xll(x — 1) ; accordingly when x is a positive integer U[x) — x ! , 

Since k'~ is less than unity 

1 



A 



oo 1 3 2^1 

2 1A111^„ K ^sin^ 



X 







r 



and since 



C\rr 1 3 2r-l JL 3 2n-l 

cos 2 * sin 2 '* 0d0 = rr l^l^lH^?^!'"^ 



*|ir 






?& 4- r I 



t% 



A 



n 



X 



/JL 3 2J--1V-2 /; 

\2 * 2 • • • — 2~~/ «- 



•* Mr. HoBSON kindly gave me^this proof when I had shown him the series which I believed to 
hold true. 
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or, with the usual notation for hypergeometric series, 



->nl 



This series is of no service, since it proceeds by powers of k*, which in our case is 
nearly unity. It is required then to transform the series into one proceeding by- 
powers of k 2 , 

It is known that, if k z + k 2 = 1, 

i< (a, 6, c, k *)= — — — — - — - /< (a, 6, l + a+6 — c, *r) 

x II (c— a— 1)11 (c — b — 1) \ - - • / 



-fr 



+/c 2(c ~ ft ~ &) — J77 — '^twtt^ — Fie — a, e—h, c — a — 6 + 1, /c 2 )." 

11(^ — 1)11 (6 — 1) x " ; 

If we apply this theorem with a = 6 = J, c = ?i+l, the first jF becomes 
7^(1 ^ i — n ? /^ whose nth and all subsequent terms involve zero factors in the 
denominators. Also the coefficient of the second F involves II ( — n — 1), which 
has an infinite factor. Hence the formula leads to an indeterminate result. Let us 
therefore put c = n + 1 + e, and proceed to the limit when e = 0. 

We have then 



7i—-€„ fC^ 



+ k^ ; j=-r— yr^ jP (n+f+e, n +■§ + €, n+l + e, k z ) >. 

Now II ( c )= 1 + € n ; (0), n ( — i+ € )=n (—J) ( i +e ™— ). 

\ ai v — j)/ 

Therefore, when e is very small, 

II (u — l-fe) / 1 1 \ 

- nT«-i) = " re V «-T + ^2 + • • • +*+ 1 )+ c J1 (°) 

= 1 + € \ II (0) + J§ -■ — — J 

HJ>^2 + e)__. s / 1 1 1 \ ir(-|) 

U(n— J) M - jC U»-1 + 2»-3^ •"• +3-1-1 ) + € n (-i) 

/T~X f ( ~»\\ n 1 \ 

ii^ / A" ~%' \ q x* \ 

. K XX \ 9/ <£jt/ ""*"' X / 

* I have to thank Mr. Hobson for giving me this formula, and for showing me the procedure wherehv 
it can be made effective. 
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Hence for the coefficient of the first series we have 



U(n-l + €)U(n + e) U(n-l)Tl(n) 



[H(«-He)P 



n(«-i) 



/ n i I r\ f ( i\ n i 






But 



n (-£)=*!* 



n'(o)-f^=log.4* 



[n(»-£)] 8 =ir ; 



^/6 i 



22» 7;, ! / ' 



II (n) IT (n — l) = n ! u — 1 ! 



Therefore 



7T 



2n! II(tt — l + e)II(tt + e) o9 r/,!^-l! 



2^(n\f [n(^~-i+6)] 2 



2^ 



1 



1 



n 



2n\ 



l + € 2 log 4— - + 2T- — 42^ — ? 



/ -» 



This is true from n=oo to 1, but in the case of n = we have 



so that in that case 



n(-i+e)=-n(e)=-+n'(o), 

II (-1 + e) 11(e) 1 
7T — ~ LJL — — = --4-2 loo- 4 



Now consider the coefficient of the second series. 



We have 
and since 



ic 2n+2< = ic 2n (l + 2clog tf ic), 



n (— x) n (a? — 1)=— 



7T 



sin 7r.^ 



II ( — n— 1 — e)U (n+e) 



7T 



( — \ n+1 



sin (% + 1 + e) 7r e 



7T. 



n(-i);r=:7y. 



Therefore the coeffi 



2n I K^ n 



icient of the second series is ( — ) n+1 -^r~r L ^ — (1+ 2e log k), and 



D = 



2?»n\ Ti^-ll 
2nl 



l + € (2log 4-^ + 2^^-4^^ 



n 



F(\, h l—n — € } K 2 ) 



+ (~Y +1 2^ ^ i (!+ 2€ lo g«) *>+£+«> ^+i + e ? w+l + «, k 3 ). 



The case of n=Q is an exception, for the coefficient of the first F has the part 



1 



inside [ ] replaced by - (1 + 2e log 4). 



* 



Proved by differentiating the known formula II (#«-!) II(ic-J) = II "(2a;- 1) . ■-■— i and putting 



r - 1 
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It remains to consider these two F series. 



r ^, s , i e,K '~^2 2 '(l-*i-€)(2-?i- e )...(r-»- e ).7! 



•vl—l 



— -^ \ / O o ; . 



[1.3...(2r-l)] 3 ^' 



(?i — l+€)('/i— 2 + e) ... (n—r+f)/' ! 



When rO 



,/ V y _ i[1.3...(2r-l)]»^ 

"^ ' ^22'(w / -l+ e )...(l+6) € (i«-e)...(>'-7i--€).r! 



1 



«- i 



(?i — 1 -f- e) (7?, — 2 + e) . 



1 / '"~ A 1 

(7i — 7* + e) "" (n—l)(n.~-2). t .(n — r)\" "^ t I 



«—)• 



When r>n, put r — n-fs, and 



1 



1 



1 



(n — 1 + e) (71 — 2 + e) . . . (1 -f e) e (1 — e) . , , (,9 — e) n — l\$l * e 



11 1 "1 € 

1 — ejj - ~+ e > - +- - 

1 1 



Also when r~n-\-s 



K" 



[1.3...(2r-l)] 3 _ _{2iilf 



92/ 



* 1 • 



/< 



2>^ 



'(2» + 1) (2»+ S) . . . (2ra + 2,9- 1)]3 

1 — — — • — — — — — — — t/ , ' ; - t 



2**(n + l)(n + 2)...(n + $) 



_L IIUS 



^ \~2> 2' ■*• ^ ^J ^ / 



r =3t-l 



y / y [1 . o ... (2?'— 1)J" / 



•n-l 



1 






Br 



n—r 



,„ (2ft !) 3 
! 4 »(w!)' 



$ = e© 



s=0 



[ (2-/1 + 1)... (2ti+ 2 g — 1)P 1/ 
2 3 *(u+l)... (w+s)s! e 



91 



+(-)»;^&^^ k ^ ;(i-^:H: + ^;-)^ 



1 . e 

t n 



* 1 



1 



It follows that we may write the first term of I) as follows :— 



(_)« — 



2ti! 



fC 



in 



2 2 "(« !)- 



^'/i : „ Z-' (71— I) ... (71 — 7*). 7' ! 



27i ! 



00 



^ V } 2^ (n If A 







[(2n + l)...(2«, + 2* -l)] 3 
2 28 <> + l)...(w + s).s! 



1 



1 



% 



1 






+ J5*y — 4^--— +2 log 4 



2« 
K. « 



The first of these terms becomes infinite when e = 0. 
Turning to the second F we have 



p, , ! , . j . , , . ox * [27i,+ l + 2e)(2»+3 + 2 6 )...(2M + 2s-l + 2 6 )l 3 , 







2 9i> (71 + 1 + e) (71 + 2 + e) . . . (?i + s + e) . s ! 



GO 



^ 



[(2»+l)(2n + 3) ... (2vi + 2s-lVl- f. 



'/> + 



1 



71 4- 5 



2 2 *(^ + l)(7i. + 2) ... (w + .9).s! 






K 



2f 
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Thus the second term of D is 

(-) ;^f-*>+*,n+i,n+i,K») 



9, 



2»(^|)S 



27i ! 



CO 



\ v 



[(2^+l)...(2^+2g-l)] 3 
2 3 *(^ + l)...(^ + ^>! 



• '/I + S -J 

4-s 



S? + 8 



1 



*" 2^-1 



^ -+2log/c 



AT . 



The first term of this becomes infinite when e = 0, but it is equal and opposite to 
the infinite term in the first part of D, and they annihilate one another. 



]j ._ _ g v_ ._ — v. T / : K ~, 



2?i\ 



o 



2* r (n~l)...(n-r)r\ 



2n\ 



•X 



/ v*— — -- k 2 ' 1 y l>~— 

V 7 2 2b (7?-!) 3 " -^'^ 



|-(2^ + l)...(2^ + 2g--l)] a 



4- 



»,+ 



1 



il+, 



2 log ~+X y +X;r - 4 2 j t 



1 



— i 



2*-l 



AC' . 



On examining the case of ft = we find that this formula also embraces it, 



provided we interpret X as nieaning zero. 

i 

The coefficient in the last term admits of some simplification, for 



n + s 



y -i l. y : *. — 4 V — — 

^t j. i <** j ^ Of 1 



n + $ 



1 



t-il 



2y . — 
A *(2*-i.) • ^ /+* 



**" i 



/-I 



1 



We thus conclude that I) or 



/"» .w 



3OS 2w 



70- 



2^^!^ — 1! 



2w! 



P 



2 2 (-w— 1)1 ! 



2 i 



J." . <->" 



.'IT 



'* %l * 



L / _ \ _1_ ^2* 



S*(?? — 1)0^ — 2)2! 



4 ''' 1 '' 1 

2 loo; —^ -; —2 2;;~ 



/c 4 — . . . to n terms 






+ 22 ^-| T)1 ,(2 log - k -x 7—2^7 



1 



(2/ - I ) ) 



K 






. . (84). 



The second integral E may he found as follows : — 

E„= (cos 2 " &dd = J cos 2 " [k"+(1 — k") cos 2 



cW 



K"±,) )l mm j m ( I "~~" K~ ) XJj^.j * . 



« » # 



* . I O «./ / • 



YOL, OXO VII. — A, 



3 Z 
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From this 1 find E or 

" «„ a a m &" +1 (n!f , &"~ l *> l n* 1 l '£ / v 3 • 3 " - (2*— l) 3 (2r+l) 
cor** . Arf#=:~~^Vf + — STTTT— X ( -)' ^7 



2??- + 1 ! 



2?i + 1 



o 



(n — 1) (w~ 2) . . . (n — r) . r + .1 ! 



K 



?r+2 



(-)» 2*i ! ? 1, (2w + 1) (2n + Sf... (2tt + 2a-3) s (2h, + 2s-1) 

_J~_ ' «/*''-•• V . — - — 

r 2 2,,+1 n \n+l\ ^ 



2 2 *- 3 (n + 2) (vi + 3) . . . (w + s) . * - 1 ! 



X 



H + .9 i S — I -1 'H + ,^—1 



log • +^ 7 + ^ . 



4 V - 



2/-1 2n + 2s — l 



K 



o .-. 



i»--l 



This is applicable also to the case of n = 0, provided that X * m interpreted as «ero. 



Id the particular case in hand I find, however, that it is shorter not to use this 
general formula, but to carry out the transformation (85) in the particular cases 
where the result is needed. 



§ 20. Reduction of preceding integrals; disappearance of logarithmic terms. 



In the application of the integrals of the last section, we are to put tc 

and only to develop as far as fi\ 

Then to the proposed order k 2 = 2)8(1 — /3), k 4, = 4/3 2 , 



1 -£ 

iTjs' 



Also 



4 8 8 

2 log -—log -5+ log (l+i8)=log 5+j8— 1^ - 



It will now facilitate future developments to adopt an abridged notation, 

then 

2* n n\n-l\ 



I writ 



e 



/ ( % )= 



9-n f 



*mA 



■and observe that t /*(?? + 1) 



2n 



2n + 1 



? /(n),and/(l) = 2, /(2) 



4 



Since k 2 is of the first order in /3, only the first series in the D integral (84) enters 
when n is greater than 2, In that ease 



D=/(n) 

=/ (h) 



/3~/? : 



_j__ 



!>/3 :! 



2(w— 1) ' 8(n-l)(w-2)_ 



/8 






(4% + 1) /3 2, 



2(7i — l) ' 8(?i-l)(tt,-2) 



# * ti 



(80). 



This result may be obtained very shortly without reference to the general formula ; 
for when n is greater than 2 



D = 



{1 + 

\1~0 

/t+y3\- 



\ 1- r>l, 



cos 3 '* eZ$ 

4 . (co^?T2^2/3^ ? 






L 



/3 + /3 3 3/3 3 ' 

'oos 3 #^2cos 4 



d0. 
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The integral of an odd power of cos 6 is easily determined, and it will be found that 
the result (86) is obtained. It is, however, clear that if n is not greater than 2 the 
development in powers of sec 3 6 is not legitimate. 

When n is not greater than 2 the formula (84) of the last section is necessary, and 
we find 



cos 1 6 



8 






- a 



^=^logs+/(2) 



13 



til 



W/8----/8* 



16 



'cos 2 6 

(18 . 

A 



CW: 



(P+lp)log%+f(l) 



l+^ + J> 



"> ( 



* e 



(1 +^+i^) log |-^ 



. .(87). 



) 



We have now to find the second integral E, and this may be done more easily than 
by reference to the general formula of the last section. 
We have 



ox fcos~' H: * 6 



fcos 3 '* 6 



E=fcos*" ^^=(1-^)]^— - tW+K 3 |^V- d$, 



A 






A 






It will be observed that even when n is 2 the general formula (8.6) gives the 
D integral as far as the first power of )3. Hence in finding E we may use that 
general formula except when n = 0, 1. 



Then since f(n) 



2n + 1 



%i " 



/* (n+1), when n is greater than I , 



E=/(n+l) 



ii-W+sp) (i -£+&£>£ )+-—? d 



2'H 



8w(w — 1) 



^i 1 ' a («~i)y J' 



1. + : 



(4ft- 1) 



2m 8«- (m— 1) 



P* 



#•*••«••«* 



* «. « 



(88). 



8 



8 



But when n = 1, 
E = (l-2j8 + 2^)^^Mog^-2i8(l^)^(l + f/5)log| 

+/(2) [(1-2^+2^) (l-^-W 3 ) +3^(1-^(1+1^); 



8 



=-^io g -|+/(2)[i+^+H^: 



. ( 89 )< 



8 



8 



And when w = 0, 

E= -( I -2/8+2/8^/8(1 +|/8) log-- +2/5(1 -/5) (1 +|/5) log ^ 

+/'(!) [(1 - 2/8+2/8*) (1 + I/5+H/5 3 ) - 2/5 (1 -/3) . rt/Vl 



8 



/8(l-i/8)log— +/(l)[l-^+i^J . . 

3/2 



* * * 



* * * 



(90), 
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I now wish to show that, in the use to which these integrals are to he put, the 
lop-arthmic terms disappear. 

O X X 

The following Is a table of these integrals collected from. (87), (89), (90), in as far 
only as they involve logarithms :— 



lW ( i +1/3+ A/8-) log j , \±de=fi(\-tf) log-- ■ , 
tW= ~p (I + */8) log -~ , (cos- A </0= - 1;/3' ; log ~- 






*cos° # 



"cos 1 6? 

A 



cW^^/3 2 log--- . 



Then the formulae (88) for L, M, N, in so far only as is at present material, are 



L 



A 



(COS 1 #+/3/c0S 2 6)-\- l r- jO/cos* $ 









« 



/ft 



.4.1 



o 

N 
H 



y A (cos 3 #+/3<7) + y. 2 A/3</cos~ 0] d.0. 
J- (cos* +/ gfc)_|_-y ( cos * 0+ph cob* 8)+°%fih cos*0 






A 



cW. 



On using the integrals and only retaining squares of /3, we iind 



L x 

?=0Vo (I-/) log ^ . 



I 



,T = £Vu (ff-^) io 8' 



(S 



X 



)8 



H 



{0*o [-(l+&8) + A (l+^)J+^a 2 ^-A)}log- 



s 







But by definition of/ and e/ in (80) and of the an in (82), to the order ^ero of 
small quantities, 






gr=J-, /^=l+i^/+-4"iS, a =-l» a 3 =— ^+1) 



&/' 



Thus the logarithmic terms entirely disappear, and henceforth may be dropped, 
Thus, as far as material, we have the following table of integrals :— 
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. A " 
cos* 4 8 



i%P" 



f cos- 



^=/(.L)Ll + f/3+W- 



C 



0S 1 M ^ 



^=/(2)[l-^-H^j. 



d0=/(fi) 



1 — 



/3 






(4» + l)/3 £ 



2(ra~1.) 8(h-:I.)(h-2)_ 



* f *•' ^***^ •"-^ * 



A^=/(l)[/l-i/3+^J. 

'cos 2 0At/0=/(2)[l +&8 + W*]. 



cos' ; "#A^0==/'(%+l) 



X |" t 



/3 (4m-1)/3 ; 



2k 8»(u- 1) 



, u> 1 



> 



/ * 



4 » 



/ 



(9.1). 



Before using these for the determination of L, M, N, it is well to obtain one other 
result. 

We have seen in (82) that 






cos 2 *' 0. 



m 



JL Ilol t/lA/A t> 



+ 1 
-1. 



(P)^= ^±i}i ( -y+-/(r+ 1) _ ~~vr- 

v 7 ^ *— si , v y t7 v ' 2 3 *V + s ! (r !) 3 7* 



6* i 



2 i 4- s ! 
I Jut this integral is equal to -.-~-~r t — ;; therefore 



M~~)' + t/ V r + ] ) 2^7+7T^r^7! = 2i + l * 
Putting #=1 and 0, and comparing with the values of a 2n y 2 ,— a * n (82), we have 







o&./O' + l) 



V> 



-^ 



iy,,_,/(r+t) 



2i + l ' 



S. 



2^4-1 i-1! 



9 f t * t * \ J * 



§ 21. Integrals of the squares of harmonics token s = 1 and s = 0. 



In (83) we have 



I 1 i 



; ri-jr 



■ -'(TT 



^(gOH 2 '^ 2 $ + /3fGOS 2 '' 0) fW. 



2r-f 1 . 



Therefore, noting that f(r)= ^77 : /( r H~l)? ailc ^ using ^ ie integrals (9 L). 



% 



+ro/(^)[i~^-^'+l/S(i+f/3)/j. 
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Substituting for f (which the reader must not confuse with the functional /' in use 



ft 



here) its value (80), the terra of order zero is ty 2r ,^ 2 f(r~{~ 1). By (92) this is aqua 



O 



to 



t> ~j~* -I, * 



2-i+l i~V. 



? • 



The term of the first order in $ is 



1 i)(2f4~ I) 



zr 



which may be reduced to the form 



/3ty 2 ^ 2 f(r+l)+^ty 2 ^^\r) 3 and is equal to -y— ^~T7+^2y 2 r-2/('' , > 



The term of the second order in (3 is 



^ty 2r ~2f(r+l) 



4r + 5 3(2'/- + 1)/ 



Vo/(2) 



+&/+£ 



61 I 9 /Si 1 ; l IV 
I 6 S 4 Vf ~T T4J "T 6 A 



This may be reduced to the form 



/3 2 






2 i + I f 
of which the first term is -|/3 3 * £77-7 7~~r» * 



m 



fherefore 



.Li 



V4. ~|~ j <£ — . J 



2 * + l! 



p^fl+zS+^g^^+i/SCl+i/SO'+^MSy^/W-iiSVo- 



N 



ow 



y»- 3 /('-)=(-)" +l • 



6+ X i 



i + rj 



It r + H-r! i— >•! 



f 5 



H 



nd 



4(-) r+1 * + *•!(*+!>• (*+l)(i+2)i(*-l) 

tr+T!»'!t-r! 1!2 



2! 2.3 



+" « » 



It is known that 



F(ci> b> e, 1) 



1 — F(?H~1 -, — a, 2, l) f 



c-l)n(c-^-6-l)* 



n(«-/<~i)n(c-.&-:i) 



* # » 



(93) 



Then since II ( — *) contains an infinite factor 



F(i+'l y —i) 2, 1) 



U (1)11(0) 



therefore 



y 



fy*r-*f( r ) 



n(-i)n(i + i) 

i + 1 ! 



0. 



« * « 



« * * 



. » ( y o )> 



* 



I have again to thank Mr. Hobsox for this' formula, which is due to Gauss. 
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ii\.lS 



O 



Hence 



L 2 i + 1 ! 'i + 1 ! 

F = 27+1 f-fi ( 1 + ^ + a^ 2 ) + i z 8 £Zjj ( 1 + t£)- 



Introducing the value of F as defined in (79), we have L or 



fai&wydv 



i + 1 ! 

2i + l ' i-T~! 



:i-P+alo^(/~12,/+68)] 



+ ff M^j[^-^} 



We have in (74) obtained p^* 1 ) 3 ^^* 1 ) 3 — (P* 1 ) 2 ]^? and if it be added to our last 



result we see that the term which does not involve the factor 1/(2^+1) is annihilated, 

and 

2ttM. i + ll 



fp(^ 1 ^ 1 F^=|^ |r^{ [l+i^(y-2)+^ T /3^(^-26/+48)] . . .(94). 



Now from (37) the square of the factor for converting & 1 into S 1 is 



"1 



' (oXIJ./ 



2. 



l+^+TfeW-8). 



Therefore 



^(f^S.T^: 



2ttM f-hl! 



2i-fl i— 1 1 
These are two of the required integrals. 



1+1)8 0'+ 2 )+iT8Ti8* (/+IQ/-96) 



» « « \ * / * 



Next Ave have from (83) 



^ = k L,_ 2 (cos 2r 6+ fa cos 2 "" 3 0) Ad0. 

VT 1 J 



2v-+l 



Noting as before that/(r) = --— — : /(r + t), find using the integrals (91), 



2? 



G 



+7o/( 2 ) 



L ^2r 8(r~l)r^ 2r py \/ r 2(f-l) / 



i+P+H/3 2 +t/%(W£) 



Substituting for # its value from (80), I find the term of order zero to be 
^72.-2 / V r i" i /' 01 2i + l ' i — 1 ! * 



The term of the first order is j8 2ya,.-j/(^+l) 



*1 2r+r 
2r 4r 



+/8yo/(2) (-Hi). 
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This may be reduced to the fovm — ft ty 2r „ 2 f(r +l) + %j3ty 2r _. 2 f (r) ; which by (92) 
and (93) becomes — ■— — ~ — i- t -\-*B -.- 



+ 1 i+l 
The term of the second order is 



1 ! 



/&i 7zr _ 2 f(r+l) 



8r(r-iy 4r \2(r- V)^ hJ ^ 2 



+ ^y /(2) ( H+^V,/ 



This is reducible to 



I/3 2 %,- 2 /(r+ O+Tff^ 2 0-4) % 2i ._ 2 /•(>•)+! £Vo : which becomes 

1 1 ' 



w 



riierefore 



2 m i -/3+i^)+ ^ J ; [p+i^ o- 1 )]. 



G 2i + 1 i-1 



?:-.1 ! 



Introducing for G its value (79), we find M or 



r** 



* — ij'TT 



h^'^JG^—^ 1 ''' 



2ttM i -f- 1 ! 

27+7 7- Tl 



l-^+ 2 l#(f- 12/+4) 



-f- 7rM. w ; 



i 4- 1 ! 



1 ! 



3 



|0+ij8*(2./+7) 



But in (74) we have ^(Ci 1 ) 2 T(P; 1 ) i! -(.P l - 1 a/^™ 8 ^^)" 



cos 3 6 



da. If this be added to 



the result just found the term which has not ]/(2/+l) as a factor is annihilated, and 



V7T 



,,.j Iff ^J t. 



*> --'t.TT 



2ttM i + 1 ! 
+ 1" 7 r ~ll' 



I + ^(3/+l0) + 3i4^(29y 2 + I34/ + -!84) . (95). 



Now from (37) the scpare of the factor for converting (£' into C 1 i« 



D, 



i (cos) 



1-J/8+AW-8). 



m 



rherefore 



iff 



2ttM i + 1! 



vTT 



p(P < 1 O < 7rf < r==r fT ^^[l+f/80-+2) + l dl^(29/+74y+48)] . . (95). 



These last two complete the solution for s=l. 
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Next we Lave from (83) 



- = ]£ 

i 1 o 



«2/ 






~{<xxi i '-+*O+ph<M?'-6)d0. 



Proceeding as before. 



N 
H 



= 2 «•.,/( r+ 1) 



2r ^ 8r (r -1) ' ^ 2r ' l -0'-" 1 )/. 



+«*/(2) [1 ~-i/3-f*/8»+&M (1 +P)J 



Substituting for h its value (80), we find that the term of order zero is ^ia 2r f(r-\- 1), 



2 
and by (92) this is equal to rr 



2i -f 1 
Die term of the first order is 



m 



/32a 2 ,/'(r+l) 

2 

which may be written in the form 



1 , 2r + r 

2r 2r 



+!/3« 2 +3/3a 



0! 



2yS 



/E2a 2 r/(^+l)+£a , and is equal to rf— :+£. 

2'^ ~r I 



o 



The term of the second order is 



/3^«,,/(r+l) 

2 



4r + 5 



- 8r(r — l) 4r(r — 1) 






8r 



+^%(U-1%) + J ^- 



03 



which is equal to 



N 



ow 



2 / v — v 2 

«r|^r i/(r+1) ~ f/(r)+i -i /(r) - 



Hence the term may be writter 



i-L 



r+1 ,, 



2 1 ' 2 



But 



And 



u 2r f(r) = (-y 



1 i -f r I 



(r !) V i — vX ' 



r+1 



aa . +2 /(r+l)=-(-^X 



1 



(r!) 2 (r + l) 3 7 



~[i(2 + l)- r ( r +l)] 7 



'i-fr ! 



£ — r I 



(->• 






i + r ! 



VOL. OXOVII.— A, 



(r!>V* (r+1) 3 ' (r!)2(r + l)Ji-^! 

4 A, 
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In the preceding formula the sum of this last function had limits i—1 to I, but as 
we j iow see that it vanishes when r=i, the upper limit may be changed to i. 
It follows that the terms of the second order are 

p4a,./(r+l)-|^/(2)-l^a /(l)- ' 



+i/3"S 



(r I) 3 ?- ' <V + 1) 3 ' 0'!/(7-+l) 



3 , (-)" 



* + •?'! , , m •!*( — Y i + rl 



The term in % in this expression will be found to be — fa 2 . That in a will be 
found to be Vce , Then since a 3 =— -^j, a =l, these terms are together i/8 2 (5/+l"0* 
The whole may then be written 



tf^f(r+i)+mJ+i)*J~]\ < J—i+i^ATTT^ i ± ':i+^(5/'+ u). 



N 



o\v 






i r I r-f-1 ! '£—r ! 



(i +!)£(•£+ 1) (i+ 2yi(/t— I) 



70'+l!>*i-H 



11 2 






i' 7 ( i + 1 , — i 5 2 , 1 ) — 1 = — 1 



v 



£ (_)r i + r \ 



i(r + lf i 



n 



*(* + !) , (i + l)( i + 2)i(i -l) 

YlJl + 3 131 



1 



9 f 1 9 "^ *J | 1 9 9 *» * * 

«W « JL ft W t^ « JL * mLi ft *L/ 



J F(i 9 — i — 1 , I , I ) — 1 + i (* + -I ) 



1 1 

:iJW; — {— 1 1 1)+ T ^1 



i 
j 



l. 



The last result follows from the fact that in accordance with (93). the sum of the 
hypergeometric series has an infinite factor in the denominator, and vanishes. 



*3 



Then since by (92) t a 2r /(r+ 1)= 9 ■ ~r . r, the terms of the second order are found to be 



2 



2M * 9,; i 1 i&PJ i" T P s 



Hence, collecting terms. 



H 



(l+iS+^^+iS+i^o'+ii). 



Substituting foi' 11 its value (7D), we have N ov 



4ttM , 



i-tt ^ ~r i- 
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But we have already found m (77) the value of | p (C*) a [(^i) a — (Pj) 2 ] dcr, and on 

adding it to the last result the term independent of l/(2i+l) disappears, and we 
have 



Vi-TT 



p(^C:) 3 ^==|^[I+Py + A^(^-10./)]. . . . (96). 



— W7J- 



2* + l 



The square of the factor whereby (£,- is converted into C; was found in (38), 
namely, 

;s=i+iW-8). 



Hence 



"v/TT 



47rM 






(96). 



These are the last of the required integrals. 



§ 22. Table of Integrals of squares of harmonies* 

In this section the results obtained in (71), (72), (94), (95), and (96) are 
collected. 



* After having completed the evaluation of all these integrals, I found that they may be evaluated 
very shortly by means of the factors <£ and E of (48), § 10. 

I find that for all values of s (writing the eight forms in a single formula), 



v f (»•)* x («i') s or (C*')3 lW __ 4*rM , , yi 



«<* 



E-* 



. x < con 



+ . r (Cr s or Ci* or §fri* or Sr v ) 2 1 

St. Dart Of V ~-~ y-— yv — ^-.__2ii_/ J.. 

(1-/3 cos 2<f>)* J 



I leave the reader to verify that this is so. 

Unfortunately I have hitherto been unable to prove the truth of this except by the laborious method 
in the text. I do not therefore know whether the result remains true for higher degrees of approxima- 
tion, although I suspect it does so. If it should be true, it would be very easy to compute the integrals 
when higher powers of /3 are included. 

It may be worth mentioning that the variables are separable in the integrals. Thus, when $3;* £i* 
denotes any one of the eight forms, 



J \,n,3. S(r ., 12 , /(r i V'WYde 










(i~i8cos 2<i>y (&i*y <i<i> 



vp 7j 



l+(i 



A "' Jo (1-/3 cos 2</>)' ''' 



— :,tc 



The c/> integrals present no difficulty, but with regard to the others we are met by the impossibility of 
expanding in powers of sec 2 for the lower orders. It would be a great step in the right direction, if it 



could be proved that all the terms which do not involve the factor 

4 A 2 



% +T 



necessarily vanish, 



548 PEOFESSOE G. H. DAEWIN ON ELLIPSOIDAL HARMONIC ANALYSIS. 



1 + flV 



It may be well to remind the reader that M=F^(^ 2 — 1)' (z; 2 — ~ — - ) 



t 



i(i + l ) 

(i — l)i (i + 1) (i -f 2) 



y=i(i-f-i ), 



First when s > 2. 
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/* 



Secondly, when $=2, 2=-^({+l). 



Type EEC f p(ftfy)G,W^=l^£Ti ^ -|/J2+ ¥ |«/? 2 (952 2 --982+72)}. 
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{1 + |)8 (S + 2)+^oj8 2 (29S 3 +1062+13fi)}. 



Thirdly 5 when 5=1, 2 is infinite and we must use ji = i(i+l). 



Type OOB \p(%W *M))*<fcr=|g£ gf ! ■! 1 +P (,/- 2 ) +T L^(/-26/+48)} t 
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Type EOS [p(»«V)S,W) 9 ^=|^^-| {l+^(j+2)+ s ^l3%f + Wj-96)}. 



Type OOC j p(P/(fi) € l 1 (^>) f g 



I 



(T 



2-rrM i+ll 



ih JZT\ { l + »& W+ 10)+ *MW+ 134?'+384) | 
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Lastly, when s=0 ; 2= — » (^+1 )= — ,/• There are only two types — 
Type EEC fp(^(M)CiW)^=|f 1 [1+Py + (^/3M7/-10i)]. 

TypeOEC fp(ft(/i)C< W)**r=^ [l + py+-i-^(7/-2/-24]. 



XxY.JlvX xX.JL. 



Summary. 



The symmetrical form in which Lame presented the three functions whose product 
is a solid ellipsoidal harmonic is such as to render purely analytical investigations 
both elegant and convenient. But it seemed to me that facility for computation 
might be gained by the surrender of symmetry, and I have acted on this idea in the 
preceding paper. 

Spheroidal analysis has been successfully employed where the ellipsoid is one of 
revolution, and it therefore seemed advisable to make that method the point of 
departure for the treatment of ellipsoids with three unequal axes. In spheroidal 
harmonics we start with a fundamental prolate ellipsoid of revolution, with imaginary 
semi-axes k*J — 1, Uy/— 1, 0. The position of a point is then defined by three co-ordi- 
nates ; the first of these, v, is such that its reciprocal is the eccentricity of a meridional 
section of an ellipsoid confocal with the fundamental ellipsoid and passing through 
the point. Since that eccentricity diminishes as we recede from the origin, v plays 
the part of a reciprocal to the radius vector. The second co-ordinate, fi, is the cosine 
of the auxiliary angle in the meridional ellipse measured from, the axis of symmetry. 
It therefore plays the part of sine of latitude. The third co-ordinate is simply the 
longitude <f>. The three co-ordinates may then be described as the radial, latitudinal, 
and longitudinal co-ordinates. The parameter k defines the absolute scale on which 
the figure is drawn. 

It is equally possible to start with a fundamental oblate ellipsoid with real axes 
k, k, 0, We should then take the first co-ordinate, £, as such that £ 2 = — v\ All 
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that follows would then be equally applicable ; but, In order not to complicate the 
.statement by continual reference to alternative forms, I shall adhere to the first 
form as a standard. 

In this paper a closely parallel notation is adopted for the ellipsoid of three 
unequal axes. The squares of semi-axes of the fundamental ellipsoid are taken to 

be — P — -- - — P, 0, and the three co-ordinates are still i>, \i^ <j>. Although their 

geometrical meanings are now by no means so simple, they may still be described as 
radial, latitudinal, and longitudinal co-ordinates. As before, we might equally well 

start with a fundamental ellipsoid whose squares of semi-axes are F ^—,^, 0, and 

replace f 9 by £~, where £ 3 = — zA All possible ellipsoids are comprised in either of 
these types by making /3 vary from zero to infinity. But it is shown, in § 2 that, by 
a proper choice of type, all possible ellipsoids are comprised in a range of fi from zero 
to one-third. When /3 is zero we have the spheroids for which harmonic analysis 
already exists ; and when /3 = -| the ellipsoid is such that the mean axis Is 
the square root of mean square of the extreme axes. The harmonic analysis for this 
class of ellipsoid has not been yet worked out, but the method of this paper would 
render it possible to do so. We may then regard /3 as essentially less than -J, and 
may conveniently make developments in powers of j3. 

In spheroidal analysis, for space internal to an ellipsoid v 09 two of the three func- 
tions are the same P-functions that occur in spherical analysis ; one P being a 
function of v, the other of /x. The third function is a cosine or sine of a multiple of 
the longitude <j). In external space the P~function of v is replaced by a Q-function, 
being a solution of the differential equation of the second kind. 

The like is true in ellipsoidal analysis, and we have P- and Q-functions of v for 
internal and external space, a P-function of /x, and a cosine- or sine-function of c/>. I 
will now for a time set aside the Q-functions and consider them later. 

There are eight cases to consider (§ 4) ; these are determined by the evenness or 
oddness of the degree i and of the order s of the harmonic, and by the alternative of 
whether they correspond with a cosine- or sine-function of <£. I indicate these eight 
types by the initials E, 0, C, or S— for example, EOS means the type in which i is 
even, ,9 is odd, and that there is association with a sine-function. 

It appears that the new P-functions fall into two forms. The first form, which I 
write :jP/, is found to be expressible in a finite series in terms of the P/ ±2I % where the 
P's are the ordinary functions of spherical analysis. The terms in this series are 
arranged in powers of /3, so that the coefficient of P/ ±2/J has /3 & as part of its 

/ 2 1 

coefficient. The second form, which I write Pf, is such that A/v — iT/s P* (*>) ov 

a/ t+f~~ P* (fO ' m expressible by a series of the same kind as that for IB .\ 
Amongst the eight types four involve \$ -functions and four P-functions ; and if for 
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given 6* a p /-function is associated with a cosine-function, the corresponding P/ is 
associated with a sine-function, and vice versa. 

Lastly ? a ^-function of v is always associated with a |9 -function of /x ; and the like 
is true of the Ps. 

Again, the cosine- and sine-functions fall into two forms. In the first form s and % 
are either both even or both odd, and the function, which I write <£/ or %\, is 
expressed by a series of terms consisting of a coefficient multiplied by /3 & cos or 
sin (s ± 2k)<f). In the second form s and i differ as to evenness and oddness, and the 
function, written 0/ or S/ ? is expressed by a similar series multiplied by 
(1 ~ /3 cos 2(j))K 

The combination of the two forms of P-function with the four forms of cosine- and 
sine-function gives the eight types of solid harmonic. 

Corresj)onding to the two forms of P-function there are two forms of Q~fuuction, 

/ v* — 1 

such that <&* and Q/ a/ -— j+p are expansible in a series of ordinary Q-functious ; 

but whereas the series for |J/ and P/ are terminable, because P/ vanishes when s is 
greater than i 9 this is not the case with the Q-functlons. In fact the series for the 
Q-functions begins with Q f - or Q/, and the order of the Q J s increases by two at a time 
up to s when we have the principal or central term ; it then goes on increasing up to 
s = i or i — 1, and on to infinity. 

In spherical and s]3heroidal analysis the differential equation satisfied by P- 
involves the integer s, whereby the order is specified. So here also the differential 
equations, satisfied by ijjj/ or P/ and by (£/, <$/, or O/, S", involve a constant, but it 
is no longer an integer. It seemed convenient to assume s z — /3<r as the form for 
this constant, where s is the known integer specifying the order of harmonic, and 
a remains to be determined from the differential equations. 

When the assumed forms for the P-function and for the cosine- and sine-functions 
are substituted in the differential equations, it is found (§6) that, in order to satisfy 
the equations, ficr must be equal to the difference between two finite-continued 
fractions, each of which involves /3cr. We thus have an equation, for /3<x, and the 
required root is that which vanishes when /3 vanishes. 

For the harmonics of degrees 0, 1, 2, 3, and for all orders, a may be found 
rigorously in algebraic form, but for higher degrees the equation can only be solved 
approximately, unless /3 should have a definite numerical value. 

When ficr has been determined, either rigorously or approximately, the successive 
coefficients of the series are determinable in such a way that the ratio of each 
coefficient to the preceding one is expressed by a continued fraction, which is, in fact, 
a portion of one of the two fractions involved in the equation for /3<x. 

Throughout the rest of the paper the greater part of the work is carried out 
with approximate forms, and, although it would be easy to attain to greater 
accuracy, I have thought it sufficient, in the first instance, to stop at fi'\ With this 
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limitation the coefficients of the series assume simple forms (§ 8), and we have thus 
definite, if approximate, expressions for all the functions which can occur in 
ellipsoidal analysis. 

In rigorous expressions, |J* and P* are essentially different from one another, but 
in approximate forms, when $ is greater than a certain integer dependent on the 
degree of approximation, the two are the same thing in different shapes, except 
as to a constant factor, I have, therefore, in § 9 determined up to squares of ft the 
factors whereby P- is convertible into P-*, and 0/ or Sf into C* or &?. With the 
degree of approximation adopted there is no factor for converting the Ps when 
s = 3, 2, 1. Similarly, down to s = 3 inclusive, the same factor serves for converting 
0/ into (£,* and S/ 5 into %■, But for s = 2, 1, one form is needed for changing 
C into C, and another for changing S into jg. It may be well to note that there 
is no sine -function when s is zero. 

The use of these factors does much to facilitate the laborious reductions involved 
in the whole investigation. 

It is Avell known that the Q-functions are expressible in terms of the P-funetions 
by means of a definite integral. Hence $1- and Q- must have a second form, which 
can only differ from the other by a constant factor. The factors connecting the two 
forms are determined in § 10. 

The second part of the paper is devoted to applications of the harmonic method. 
In § 11 the perpendicular from the centre on to the tangent plane to an ellipsoid v , 
and the area of an element of surface of the ellipsoid, are found in terms of the 
co-ordinates /x, <f>, and the constant v . 

It is easy to form a function, continuous at the surface v , which shall be a solid 
harmonic both for external and for internal space. Poisson's equation then enables 
us to determine the surface density of which this continuous function is the potential, 
and it is found to be a surface harmonic of /a, (j> multiplied by the perpendicular on to 
the tangent plane. This application of Poisson's equation involves the use of the 
Q-function in its integral form. Accordingly, if the serial form for the Q-function is 
adopted as a standard, the expression for the potential, of a layer of surface density 
involves the use of the factor for conversion between the two forms of Q-function. 

This result may obviously be employed to determine the potential of an harmonic 
deformation of a solid ellipsoid. 

The potential of the solid ellipsoid itself may be found by the consideration that it 
is externally equal to that of a focaloicl shell of the same mass. It appears that in 
order to express the equivalent surface density in surface harmonics, it is only 
necessary to express the reciprocal of the square of the perpendicular on the tangent 
plane in that form. This result is attained by expressing x 1 , y°, z* in surface 
harmonics. When this done, an application of the preceding theorem enables us to 
write down the external potential of the solid ellipsoid at once, ' In § 12 the external 
potential of the solid ellipsoid is expressed rigorously in terms of solid harmonics of 
degrees zero and 2. 
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Since x 2 y y 2 , z 2 have been found in surface harmonics, we can also write down a 
rotation-potential about any one of the three axes in the same form. 

The internal potential of a solid ellipsoid does not lend itself well to elliptic 
co-ordinates, but expressions for it are given in § 12. 

If it be desired to express any arbitrary function of /x, (f> in surface harmonics, it is 
necessary to know the integrals, over the surface of the ellipsoid, of the squares of 
the several surface harmonics, each multiplied by the perpendicular on to the tangent 
plane. The rest of the paper is devoted to the evaluation of these integrals. No 
attempt is made to carry the developments beyond /3 3 , although the methods 
employed would render it possible to do so. 

When s is greater than unity, it appears that it is legitimate to develop the 

1 
function to be integrated in powers of ■j--™~- 2 - ; and when this is done, the integration, 

although laborious, does not present any great difficulty. 

But when s is either 1 or 0, the method of development breaks down, because it 
would give rise to infinite elements in the integrals at the poles where /x 3 is unity. 
However, portions of the integrals in these cases can still be found by the former 
method of development. As to the residues which cannot be so treated, it appears 
that they depend on integrals of the forms 



* w cos 3w 8 cW r ** 



a 



— |tt (1 — a/ 3 sin 2 0) 



and cos 2 " 0(1-— k 2 sin 3 fff d$ 9 



a 77 



where *" is nearly equal to unity. 

Development of the square-roots in powers of k 2 is useless on account of the slow 
convergence, and it is required to find series which proceed by powers of /c 3 , where 

K v = 1 — /C ^ 

By a somewhat difficult investigation, in respect to which I owe my special thanks 
to Mr. Hobson, the needed series are found (§ 19). 

It appears that portions of the two integrals involve logarithms which become 
infinite when k vanishes. Since, in the application of these integrals, the vanishing 
of k implies the vanishing of /3, we appear to be met by a difficulty. It is known 
that in spheroidal analysis no such terms appear, and we may feel confident that 
they cannot really exist in ellipsoidal analysis. In § 20 it is proved that the 
logarithmic terms do as a fact disappear. The residues of the integrals in the cases 
s = 1, are thus found, and added to the previous portions to form the complete 

results. 

ThS second part of the paper ends (§22) with a list of the integrals of the squares 
of the surface harmonics for all values of .9, as far as the squares of /3. 

Finally, an appendix below contains a table of all the functions as far as i = 5, 
5=5. It is probable that for the higher values of ,9 the results would only be 
applicable when fi is very small. 
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Appendix. 



Table of the P- and Q- Functions. 
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Note that in this table P* denotes 
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If the variable is /x ? and if accordingly the factor (ir — 1)** in P/ is replaced by 
(1 — fi 2 )**, the signs of all the terms which have /3 as coefficient must be changed, 
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"Fable of the Cosine and Sine Functions, 
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OOC) 

OOS) 

OEC) 

OES) 

OOC) 

OOS) 

OEC) 
OES) 

OOC) 

OOS) 



€1 



= \Z *- W ± W) \Z 3< ^ + w 



COS r f 

. 0(/> 

sm r 



ros_ 



$ 



¥/3{i+J2 s 2^-t^{ s c ^^ 







sm 



c 4 



=m^± x m\ G Z <£+ (I s ^-^ll s 5c/> 



sm 



sm 



sm 



cos , . COS A , 

sm T ' sm r 



<Ss [sin/' lbr [sm r [sm r 



Note that in this table 



<£ = (1— cos 2<£)*. 



A table of P (i/) and Q (v) up to i = 5, s = 5 is contained in Professor Bryan's 

paper (' Proc. Camb. Phil. Soc./ vol. vi., 1888, p. 297). The functions there 

tabulated as T^ (v) and U* (v) in the notation here adopted, would be P* (y) (with the 

i — s I 
factor (^ 2 — lf s ) and (~-) s - — — Q* (v). 

The formula for Q/ (v)\ where s is greater than i, is given in § 1 above. 



